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ABSTRACT 

Thermally broadened Lya absorbers (BLAs) offer an alternative method to highly-ionized metal 
lines for tracing the warm-hot intergalactic medium (WHIM) at T > 10^ K. However, observing BLAs 
requires data of high quality and accurate continuum definition to detect the low-contrast features, 
and a good knowledge of the velocity structure to differentiate multiple blended components from a 
single broad line. Even for well-characterized absorption profiles, disentangling the thermal line width 
from the various thermal and non-thermal contributors to the observed line width is ambiguous. 
We compile a catalog of reliable BLA candidates along seven AGN sight lines from a larger set of 
Lya absorbers observed by the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space 
Telescope (HST) . We compare our measurements based on independent reduction and analysis of the 
data to those published by other research groups. We examine the detailed structure of each absorber 
and determine a reliable line width and column density. Purported BLAs are grouped into probable 
(15), possible (48) and non-BLA (56) categories. Combining the first two categories, we infer a line 
frequency (cW /dz)BLA = 18 ± 11, comparable to observed O VI absorbers, also thought to trace 
the WHIM. We discuss the overlap between BLA and O VI absorbers (20-40%) and the distribution 
of BLAs in relation to nearby galaxies (O VI detections in BLAs are found clouser to galaxies than 
O VI non-detections). We assume that the line width determined through a multi-line curve-of-growth 
(COG) is a close approximation to the thermal line width. Based on 164 measured COG H I line 
measurements, we statistically correct the observed line widths via a Monte-Carlo simulation. Gas 
temperature and neutral fraction /hi are inferred from these statistically-corrected line widths and 
lead to a distribution of total hydrogen columns. Summing the total column density over the total 
observed pathlength, we find a BLA contribution to the closure density of flsLA = 6.31q g x 10"'^ h^Q 
based on 10'* Monte-Carlo simulations of each BLA system. There are a number of critical systematic 
assumptions implicit in this calculation, and we discuss how each affects our results and those of 
previously published work. In particular, the most comparable previous study by Lehner et al. (2007) 
gave nBLA — 3.6 x 10~^ h^^ or 9.1 x 10~^ft.yQ^, depending on which assumptions were made about 
hydrogen neutral fraction. Taking our value, current O VI and BLA surveys can account for 20% of 
the baryons in the local universe while an additional ~ 29% can be accounted for in the photoionized 
Lya forest; about half of all baryons in the low-z universe are found in the IGM. Finally, we present 
new, high-S/N observations of several of the BLA candidate lines from Early Release Observations 
made by the Cosmic Origins Spectrograph on HST. 

Subject headings: cosmological parameters — cosmology: observations — intergalactic medium — 
quasars: absorption lines — ultraviolet: general 



1. INTRODUCTION 

Theoretical studies of cosmological structure forma- 
tion and the intergalactic medium (IGM) predict that 
the baryonic component should develop a hot, shocked 
phase at temperatures of T > 10^ K. Adding to these 
gravitation shocks is the energy feedback from galactic 
winds, which can produce a hot circumgalactic medium 
extending 100-200 kpc from galaxies. Understanding this 
warm-hot intergalactic medium (WHIM) is therefore cru- 
cial to understanding local large-scale structure, galaxy 
feedback, and the spread of metals. Because the IGM 
in the early universe (z = 2 — 6) was primarily warm 
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and photoionized, Lya forest absorption can be used to 
trace the large-scale structures at those times accounting 
for the majority of the baryons (Ranch 1998). However, 
the fraction of baryons in the Lya forest is observed to 
change in form and content toward lower redshifts ow- 
ing to the rapid drop of the extragalactic ionizing back- 
ground and the formation of large-scale structure (Pen- 
ton, Stocke, & ShuU 2004; Weymann et al. 2001). Cos- 
mological simulations (e.g., Dave et al. 1999; Dave et al. 
2001; Gen & Ostriker 1999; Gen & Fang 2006) predict 
that by z ~ the baryons are distributed among three 
dominant reservoirs: 10-30% in cool gas (< 10** K) in 
and around galaxies, 30% remaining in the warm, pho- 
toionized Lya absorbing gas, and the remainder in hotter 
WHIM gas. These same simulations predict that shocks 
from diffuse matter falling onto large-scale structure fil- 
aments and from galactic wind feedback into these fila- 
ments have heated a large fraction of the low-z baryons. 
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While a broad spectrum of shock velocities undoubt- 
edly produces some shocked material at cooler tempera- 
tures, the WHIM phase is defined as low-density gas at 
T = 10^ — lO'^ K. Two recent simulations estimate the 
baryon fraction residing in the WHIM at 30 — 50%, but 
with considerable uncertainty (Dave et al. 2001; Cen & 
Fang 2006) largely due to the pervasiveness and effec- 
tiveness of feedback from galaxies into the IGM. Thus, 
no accurate baryon census can be made without detailed 
knowledge of WHIM absorbers.. 

Three observational methods have been employed to 
detect and obtain a rough census of the WHIM. Ions 
with high ionization potentials (> 100 cV) are likely the 
product of coUisional ionization and thus trace hot gas. 
Absorption of highly-ionized metals can in principle be 
observed in the X-ray; in particular O VII (21.60 A) 
and O VIII (18.97 A) since they arc strong transitions 
of an abundant element with peak coUisional ioniza- 
tion equilibrium (CIE) abundances at temperatures of 
(1 — 2) X 10^ K. Cooling of diffuse gas at these tem- 
peratures is slow and a significant gas reservoir is ex- 
pected. While there have been several reported inter- 
galactic O VII and O VIII detections (Fang et al. 2002; 
Fang, Canizares, & Yao 2007; Nicastro et al. 2005) based 
on Chandra observations along AGN sight lines, none 
has been confirmed with XMM/Newton spectra (Ras- 
mussen et al. 2007; Kaastra et al. 2006) and their reality 
has been questioned (Bregman 2007; Richter, Paerels, 
& Kaastra 2008). Yao et al. (2009) used a statistical 
"stack-and-add" technique looking for O VII absorption 
associated with known O VI systems to place upper lim- 
its of iVovii ^ 10 X A^ovi- However, definitive X-ray 
WHIM detections have remained, at best, controversial. 

Less highly-ionized species such as O VI provide the 
second WHIM tracer. Until now, this has proven to 
be the most fruitful technique due to extensive use of 
the UV spectrographs on board the Hubble Space Tele- 
scope (HST) and the Far Ultraviolet Spectroscopic Ex- 
plorer (FUSE). Absorption in the Li-like doublets O VI 
AA1032,1038 and N V AA1238,1242 traces gas with peak 
CIE abundances at the lower end of the WHIM tempera- 
ture range (10'''^'' K). Several large O VI absorption-line 
surveys of the local IGM toward bright AGN have been 
pubhshed. Danforth & Shull (2008, hereafter DS08) re- 
port 83 O VI absorption systems out of a total of ~650 
Lyof absorbers along AGN sight lines surveyed B.i z < 0.4. 
Tripp et al. (2008) and Thom & Chen (2008) used some- 
what smaller numbers of sight lines, but report 51 and 
27 O VI absorbers, respectively, with detection statis- 
tics similar to DS08. Two Ne VIII AA760, 770 detections 
(Tmax ~ 10^-^^ K) have been made (Savage et al. 2005; 
Narayanan, Wakker, & Savage 2009). These may provide 
a more reliable WHIM tracer, but the current statistics 
are poor. 

There are several controversies, weaknesses, and un- 
knowns that underlie UV metal-line WHIM surveys (see 
Danforth 2009). While the number of absorbers is gen- 
erally agreed upon, their interpretation remains contro- 
versial; O VI, N V, and C IV ions are produced readily in 
coUisionally ionized gas, but they can also be photoion- 
ized by sufficiently energetic photons (114, 78, and 48 
cV, respectively). The observed hydrogen and oxygen 
column densities require very diffuse {uh < 10~^ cm~'^) 



gas over long pathlength regions (> 500 h^^ kpc) to rc- 
prociucx; the line strengths and ratios seen in Lya and 
O VI absorbers (e.g., Prochaska et al. 2004; Tripp et al. 
2008; Oppenheimer & Dave 2008). Thus, some high-ion 
absorbers may not be shock-heated WHIM material, but 
instead low density, photoionized gas. 

Regardless of the ionization state, it is important to 
note that O VI absorbers are found in locations rela- 
tively near bright (i*) galaxies (< iQQhjQ kpc Stocke 
et al. 2006; Wakker & Savage 2009) as expected for 
metal-enriched WHIM gas based upon numerical sim- 
ulations (Dave et al. 1999; Cen & Ostriker 2006; Cen 
& Fang 2006). Since some degree of metal enrichment 
is required for this WHIM diagnostic (the "metallicity 
bias"), we cannot use this method to sample the pri- 
mordial IGM far from metal production sites in galaxies 
and beyond the range of metal distribution via starburst 
winds (Stocke et al. 2006, 2007). Thus, this method can- 
not possibly sample all WHIM gas. Finally, translat- 
ing metal-ion statistics into a total baryon distribution 
requires estimating the metallicity Z and fractional ion 
abundance /ion, entailing additional assumptions. Even 
the best statistics for O VI line densities {dN jdz) may 
result in a factor-of-two uncertainty in the final baryon 
budget. 

An orthogonal route to detecting WHIM absorbers, 
and one that docs not rely upon metal enrichment, is 
to detect and measure broad Lya absorbers (BLAs) in 
high-S/N spectra. Hydrogen is so abundant in the uni- 
verse that, even with a neutral fraction of 10~^ or less, 
there can be measurable absorption in Lya lines. There 
is no "metallicity bias" and no abundance uncertainty. 
In coUisionally ionized gas, the neutral fraction is a func- 
tion of temperature and can, in principle, be directly 
measured from the line width. Line width is usually de- 
noted in terms of the corresponding doppler parameter 
h = FWHM/2\/hi2 = ^/2a where g is the Gaussian 
width. Thermal width is simply a function of gas tem- 
perature T and atomic mass A\ 

6t(T) = ^/2kT/m = ^/T/QQA km s"^ 

= (40.6kms-l)T5/^yl-l/^ (1) 

where T5 is temperature in units of 10^ K. 

HST has observed over 100 bright AGN sight lines at 
modest signal-to-noise (S/N) ratios (5 — 20) with its three 
generations of spectrometers: the Faint Object Spec- 
trograph (FOS), the Goddard High Resolution Spectro- 
graph (GHRS), and the Space Telescope Imaging Spec- 
trograph (STIS). Until the recent installation of the Cos- 
mic Origins Spectrograph, the best data for a BLA search 
comes from recent STIS echelle spectroscopy, which typ- 
ically features the highest S /N combined with a spectral 
resolution sufficient to resolve H I lines of widths less 
than ~ 30 km s~^. From detailed analyses of moder- 
ate S/N STIS echelle spectra of individual bright AGN, 
Sembach et al. (2004) and Tripp et al. (2001) noted the 
presence of a few very broad Lya lines confirmed to have 
h > 40 km s^^ using curve-of-growth (COG) analyses 
of multiple Lyman lines. Richter et al. (2004), Aracil 
et al. (2006a), Lehner et al. (2006), and Williger et al. 
(2006) presented BLA candidates based upon line-width 
measurements, as did Penton, Stocke, & Shull (2004). 
Richter et al. (2006) compiled a few of the very best 



cases of BLAs. More reeently Lehner et al. (2007, hence- 
forth L07) produced a large compilation of 7 AGN sight 
lines and ~ 100 candidate BLAs based on 6- values > 40 
km s~^ (defined variously; see below). 

However elegant, BLA surveys are not without their 
practical complications, as we describe here and in 
Section 2.1. First, identifying bona fide, thermally- 
broadened H I systems is challenging. Absorption fea- 
tures arising from trace neutral hydrogen in WHIM 
should exhibit broad profiles and low column densities, 
which are difficult to detect in spectra of finite S/N. 
Additionally, differentiating single, broad features from 
blended profiles or those with ambiguous component 
structure is central to their identification as BLAs. Be- 
cause of these ambiguities, we provide an independent 
opinion on many previously studied sight lines. Second, 
inferring a temperature from a measured line width is 
crucial for inferring the total gas column density present. 
Most previous authors recognize that the observed line 
width 6lw is, at best, an upper limit on the thermal line 
width, so that only upper limits on total baryon mass can 
be inferred. For stronger Lya absorbers, multiple Lyman 
lines allow us to measure the curve of growth (COG) 
for cool H I systems and then apply a statistical correc- 
tion to obtain better estimates on the distribution of 6t 
(Section 3.3). Finally, determining a hydrogen neutral 
fraction /hi from an inferred temperature depends on 
assumptions about collisional (thermal) ionization and 
photoionization by the mctagalactic radiation field. In 
Section 4.1, we address the relative importance of these 
two mechanisms from first principles. 

Owing to the difficulties in definitively detecting BLAs 
and the clear importance in doing so, we take a detailed 
look at broad H I systems from our own work and other 
published studies. In Section 2 we discuss our method- 
ology, including a discussion of the complications sur- 
rounding BLA surveys. Our results are presented in Sec- 
tion 3, including a cross-correlation of our BLA cata- 
log with recent large surveys of low-z IGM metal line 
absorbers (DS08) and nearby galaxies. We discuss the 
larger cosmological conclusions of this BLA census and 
its systematic uncertainties in Section 4. Our conclu- 
sions are presented in Section 5. Those interested in an 
overview of our results arc cautioned against skipping 
directly to Section 3, since understanding our selection 
criteria for BLA vs. non-BLA absorbers is critical to 
understanding our results in general. 

2. METHODOLOGY AND ABSORBER SAMPLE 

2.1. Broad Lya Lines as a WHIM Tracer 

Using thermally-broadened H I lines to trace WHIM 
gas avoids many of the biases that can plague metal- 
ion based WHIM tracers (see Danforth 2009). How- 
ever, BLA surveys are fraught with observational com- 
plications as well (e.g. Richter et al. 2004, 2006). First, 
BLA lines are hard to detect and require spectra with 
a high signal-to- noise ratio (S/N) and minimal instru- 
mental systematics. For example, at T = 10'' K, the 
hydrogen neutral fraction is approximately /hi ~ 10~^ 
(Fig. 1) and the thermal 6- value is 40 km s~^ (Eq. 1). An 
IGM absorber with a total hydrogen (H I-l-H II) column 
density A^h = 10^^ cm^^, typical of what is expected in 
the Lya forest, would still have an observable H I col- 
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Fig. 1. — The neutral fraction /hi of hydrogen drops sharply 
as a function of temperature. In collisional ionization equilibrium 
(dashed Une, Sutherland & Dopita 1993) /hi < IQ-^ in the WHIM 
temperature range (T = 10® — lO'^ K, shaded region). Richter, 
Fang, & Bryan (2006) parametrize CIE-|-photoionization simula- 
tion results to a linear relationship with log T (dotted line) that is 
~ 2 - 6 times smaller than GIF at T = (1 - 3) X 10® K. CLOUDY 
photo+thermal ionization models (Danforth et al. 2006) with pho- 
toionization parameters log U = —2 and —1 (solid lines) approxi- 
mate the simple CIE curve throughout the WHIM range (shaded), 
but deviate significantly at lower temperatures. Thermal 6-value 
for H I is shown on the top axis. 

umn of TVhi ~ 10^"^ cm~^. In the optically thin regime 
(linear COG of growth), the equivalent width Wx and 
line-center optical depth to{N,T) for a Gaussian profile 
BLA are 

W,{N) = (^) (^) = (54.5 mA) N.s (2) 

ro(Ar,r) = — ^ = (0.187) iVisT-i^ (3) 

where Ni^ is the H I column density in units of 
10^3 cm-2. Thus, an IGM absorber with N13 « 1 will 
have a rest-frame Lya equivalent width W\ = 55 mA and 
a fractional depth of ~ 20% (assuming no non-thermal 
broadening). This line would be easily observable in data 
of modest S/N. 

However the prospects get much worse at higher tem- 
peratures, as the thermal line width rises as Vt and the 
neutral fraction drops (Fig. 1). A Lya absorber with the 
same total hydrogen column as above but with a tem- 
perature of 3 X 10''' K (near the peak CIE abundance 
of O VI) would have /hi ~ 10-^ iVni ~ lO^^ cm'^, 
6 = 70 km s~^, W\ ^ 5 mA, and a fractional depth 
To ^ 1%, impossible to detect unless the spectra have 
S/N ^ 100. Turbulent and/or bulk cloud motions would 
broaden the Lya line even further and reduce the frac- 
tional line depth. Confirmation of a Lya absorber and 
definitive measurements of 6hi and A^hi using higher Ly- 
man lines and COG techniques are difficult for such lines 
since even Ly/3 is a factor ~ 7 weaker than Lya for un- 
saturated systems. 
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Fig. 2. — Comparison of line width 6lw and COG doppler parameter bcOG in 164 Lya absorbers from DS08. Top left: histogram 
of measured 6lw (dotted line) and bcOG (solid line). Bottom left: ratio 6lw/''COG covers a wide range, with a mean 1.35 and median 
1.26^q'25 (dotted lines). The top right panel shows that 6lw/''COG 's weakly correlated with 6lvv = 0.45). Lower right panel shows 
that ftLw/^COG is uncorrelated with log A^hi {R = 0.1). 



Detecting low-contrast absorption features depends 
critically on obtaining a correct model of the contin- 
uum level, as well as a presumption that the intrinsic 
continuum emanating from the AGN central engine is 
flat and featureless over these broad line widths. Since 
the UV continuum of Scyfcrts and QSOs is thought to 
arise in thermal accretion disk emission close to the cen- 
tral black hole, this requires that the accretion disk itself 
be featureless, which may or may not be a correct as- 
sumption. For this reason, the UV continua of BL Lac 
objects may provide better background sources for sen- 
sitive BLA detection. Additionally, fixed pattern noise 
and echelle-blazed gratings can themselves produce gen- 
tle undulations in the continuum, either on the detec- 
tor or in the extraction process of curved orders imaged 
on rectilinear detectors and their corresponding artifacts 
and blemishes. The possible inability to place an ac- 
curate continuum due to the above difficulties not only 
makes low-contrast line detection suspect, but also makes 
it hard to set sensible detection limits as a function of 
6- value. 

Even assuming data with sufficiently high S/N and a 
well-defined continuum, the identification of BLAs re- 
lies crucially on deconvolving the thermal line width &t 
from the total observed line width 6obs- A simple, single- 



component absorber can be modeled as a Voigt profile, 
and the observed line width is the quadrature sum of 
the thermal line width 6t, the instrumental point spread 
function (usually approximated as a Gaussian 6inst)) and 
non-thermal broadening term 6nt- The latter term en- 
compasses bulk turbulence, multiple unresolved velocity 
components, or any other condition that will broadc;n a 
line profile. Instrumental broadening is generally well 
determined and can be subtracted from the total line 
width. However, the relative proportion of the thermal 
and non-thermal contributions is typically degenerate for 
absorption in a single species. 

Some simulations (e.g., Richter, Fang, & Bryan 2006) 
suggest that turbulence contributes on average ^ 10% to 
the total line width in broad absorbers, while other sim- 
ulations suggest that it may be as much as 50% (Cen & 
Fang 2006). Observational studies have shown that the 
line width measured from a single Lya line is generally 
greater than that determined through a multi-Lyman- 
series COG (ShuU et al. 2000). For example, the so-called 
Virgo Cluster absorbers in the high S/N GHRS spectrum 
of 3C 273 show a pair of Lya components with measured 
6 = 41 km s~^ and 34 km s~^ for the 1015 km s"""^ 
and 1590 km s"-"^ IGM absorbers, respectively (Weymann 
et al. 1995). However, subsequent FUSE observations 
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(Sembach ct al. 2001) determined bcOG — 30 km s^^ 
and 16 km s^^ using three and eight higher-order Ly- 
man hnes, respectively. 

While these absorbers make good individual exam- 
ples, statistical studies show that 6-values determined 
from line-width measurements of single lines, hereafter 
denoted 6lw. systematically ovcrpredict the true line 
width. ShuU et al. (2000) found that foLw/^coG ~ 2 using 
Lya and Ly/3 measurements of 12 absorbers. Danforth 
et al. (2006) confirmed the low-^ result with a larger sam- 
ple (~ 100 absorbers). Similar conclusions were reached 
by Songaila (1997, 1998, 2001) for high-z Lya forest lines. 
We refine this relationship further using the large catalog 
of DS08. Out of ~ 650 H I absorbers in DS08 plus ad- 
ditions, we measured 164 in multiple Lyman lines, and a 
COG analysis was performed giving more accurate values 
of 6hi and A^hi- Figure 2 shows the distribution of ^Lya 
and bcoG for this sample. Note that there is no clear 
correlation of bhw/bcoG with column density and only a 
weak trend with 6lw. The median overprediction ratio 
is 6lw/^cog = l-26iQ'25 ^i^^ ^ mean of 1.35. We will 
use this result extensively in our analysis for the BLAs 
in which only Lya measurements are available. 

Another impediment to BLA identification and mea- 
surement is that Lya absorption lines often show ambigu- 
ous component structure. Several narrow components in 
a close blend can artificially broaden a line. Sometimes 
this will manifest itself as an asymmetric or obviously 
non-Gaussian line profile, but this blending can be un- 
detectable even at arbitrarily high S/N and good spec- 
tral resolution. For this reason, even COG-determined 
6-values can overprcdict the thermal line width. 

To estimate the degree of line-width overprediction 
from COG analysis, we have modeled two simple two- 
component systems where the line width and relative 
centroid separation are varied (see also Danforth et al. 
2006). In case 1 (Figure 3a), two narrow components, 
each with thermal line width 6t = 15 km s~^, are 
blended into a single absorption line in Lya and higher 
Lyman lines. The equivalent widths of each blended Ly- 
man line are fitted with a COG. For component sep- 
aration Av < 5t, the inferred 6coG overpredicts 6t by 
< 20%, but at Av > 6t: the overprediction rises to nearly 
50%. This is largely independent of the relative column 
densities of the two lines. Depending on the strength 
of the line, the resolution of the spectrograph, and the 
quality of the data, even separations of Av > fex arc not 
obvious in the line profile. Thus, a pair of narrow lines 
can easily be mistaken for a single broader system, even 
using a COG. 

In case 2 (Figure 3b), we combine a narrow system 
(^T.i = 15 km s""'^) with a broad absorber (5t,2 = 2 6t,i) 
and vary the relative line strengths and component sepa- 
rations. For components of similar column density, bcoG 
is close to the mean of the two components, and separat- 
ing the components by (1—2) 6t,i increases the composite 
&COG by ^ 20 — 50% as in case 1. When one component 
is considerably stronger than the other, it dominates the 
6coG solution. It is reassuring that the total column den- 
sity of the composite system is conserved, independent 
of component separation and relative strength (Jenkins 
1986). 

Appropriate to case 2, many observed absorbers clas- 
sified as WHIM due to their O VI absorption are mul- 



tiphase in nature (e.g., Danforth et al. 2006; Danforth 
& ShuU 2008; Tripp et al. 2008). A cool, photoionized 
component (T < 10^ K) is present at nearly the same 
velocity as a warm- hot component. We now estimate 
the broad Lya profile expected to correspond with ob- 
served WHIM gas. For a typical O VI absorber with 
iVovi « (10^^-^ cm-2)7V}°^^) at the peak CIE O VI 
abundance temperature (T « Tmax = 10^ '^^ K) the re- 
sulting BLA should have total hydrogen column density 

Ar(H) ^ ^ovi 

(/ovi)(4.90xl0-5)(g^) 

= (2.93 X 10^« cm-')Ni2r {^Y^ (4) 

and, from Eq. 1, 6hi = (68.2 km s"!) (T/Tmax)^/^- 
Here we adopt an oxygen abundance (O/H) = (4.90 x 
10"5)(Z/0.1Zq) scaled to 10% of the solar mctallicity 
(Asplund et al. 2005). Given a typical neutral frac- 
tion /hi w 1.4 X 10"*^ at T w Tmax, the resulting BLA 
will have iVm ~ 5 x 10^^ cm~^ and tq ~ 0.05, eas- 
ily overwhelmed by the narrower, stronger H I absorp- 
tion from the photoionized component (typically A'hi ~ 
lO^^-i^ cm-2, 6 = 20 - 30 km s'^. To > 1). The BLA 
signature will be apparent only in the line wings and will 
require exquisite data (S/N > 30 — 40) and a very good 
knowledge of the instrumental point spread function to 
recover. 

Multiphase systems are difficult to identify individ- 
ually, but there is statistical evidence for broad-plus- 
narrow H I systems (e.g., Tripp et al. 2001; Danforth 
& ShuU 2008; Tripp et al. 2008). DS08 compared the 
6hi distribution of 83 H I systems with O VI detec- 
tions with another 273 having clean O VI nondetections 
(A^ovi < 10^^-^ cm-2). For the O VI detections, the me- 
dian and standard deviation were bm = 31 ± 15 km s~^, 
while the O VI nondetections show 6hi = 26±13 km s~^. 
This slight difference (at a low confidence level) suggests 
that weak BLA lines might be present in the O VI sys- 
tems and broadening their overall H I profiles. However, 
it is doubtful whether a difference between the two pop- 
ulations would be apparent without the O VI detection 
"sign-posts" . There are several good, individual exam- 
ples of this observational signature (e.g., Tripp, Savage, 
& Jenkins 2000; Tripp et al. 2001; Stocke, Keeney, & 
Penton 2005; Stocke et al. 2006). 

Observing multiple Lyman lines and determining 6hi 
rigorously through a COG tends to produce more accu- 
rate line width measurements, but even here, multiple 
components can broaden a line width. Because the over- 
estimate of the true 6t is smaller using a COG than for 
Lya alone, we will adopt 5coG as our best estimator for 
absorber temperature. Examining higher-order Lyman 
lines cannot be counted on to solve the multi-component 
problem since most BLA absorbers are too weak to allow 
Ly/3 to be detected and measured in all components. Nor 
will metal ions, which show intrinsically narrower lines, 
always help, since many BLAs arc expected in regions of 
little or no metal enrichment. Thus, we must adopt a 
statistical approach to BLA verification. 

2.2. Dataset and Methodology 

With the above points in mind, we set out to determine 
what fraction of the reported BLAs are legitimately trac- 
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Fig. 3. — Blended absorption components in simulated STIS/E140M data. Two narrow components (left, 6i = 62 = 15 km s~^) 
separated by Av = &i (component centroids marked with arrows) produce a combined jHofile which is indistinguishable from a single, 
broader component even in data of high S/N and resolution {S/N = 30 per 7 km resolution element). Fitting the profile as a single line 
gives bLya = 19 km s~^ or a 27% overestimate of 61. Even a COG analysis using Lya and Ly/3 profiles gives 6cOG ~ 1-2 61. Doubling the 
line width of one component (right, 62 = 2 61 = 30 km s~^) simulates a multiphase, broad + narrow system. The asymmetric profile maizes 
it obvious that a single-component fit is not appropriate, but the component deconvolution and number of components are ambiguous. 



TABLE 1 
BLA Sample Sight Lines 



AGN 


R.A. 


Dccl. 


ZAGN 






Source 


HE 0226-4110 


02:28:15.2 


-40:57:16 


0.495 


0.389 


0.305 


Lehner et al. (2006) 


HS 0624+6907 


06:30:02.5 


+69:05:04 


0.370 


0.354 


0.282 


Aracil et al. (2006a,b) 


PG 1116+215 


11:19:08.6 


+21:19:18 


0.177 


0.167 


0.150 


Sembach et al. (2004) 


PG 1259+593 


13:01:12.9 


+59:02:07 


0.478 


0.388 


0.303 


Richter et al. (2004) 


PKS 0405-123 


04:07:48.4 


-12:11:37 


0.573 


0.387 


0.302 


Lehner et al. (2007); Williger et al. (2006) 


H 1821+643 


18:21:57.3 


+64:20:36 


0.297 


0.283 


0.236 


Sembach et al., in prep. 


PG 0953+414 


09:56:52.4 


+41:15:22 


0.234 


0.225 


0.195 


Tripp et al., in prep. 



»Most recent detailed STIS/E140M analysis; all are also included in Lehner et al. (2007); Danforth & ShuU (2008). 



ing warm-hot gas at T > 10^ K. To accomplish this, we 
have independently extracted, reduced, and scrutinized 
the STIS/E140M spectra (Table 1) used by Lehner et al. 
(2007, L07 hereafter) to search for BLAs potentially aris- 
ing in the WHIM. Complete details of our data reduction 
method are given in DS08. Briefly, STIS/E140M data 
were uniformly reduced using CalSTIS v2.19. Line- 
free continuum regions were defined interactively in lOA 
segments of the data and fitted with low-order Legendre 
polynomials. Continuum fit uncertainty was taken as 
the standard deviation of points about the mean in the 
defined continuum rc;gion. This uncertainty was added 
in quadrature with those from photon noise and line fit 
uncertainties. 

The L07 work is the most complete and comprehensive 
look at BLAs currently available. In order to determine 
the number and 6-value distribution of BLA absorbers 
in the local universe, we have adopted the same defini- 
tion as L07 and most other studies, requiring T > 10^ K 
corresponding to a hydrogen thermal 6 > 40 km s~^. 
The relationship between COG-determined line width 
and line width measured from a single line shows a me- 



dian ratio 6lw/^COG = ^■'^^-0.25 the large DS08 sam- 
ple, and the correlation of that ratio with either A^hi or 
6lw is poor (Fig. 2). Given this ratio, an absorber with 
&LW ^ 50 km s~^ (i.e., 1.26 x 40 km s~-^) would be more 
likely than not to have a 6cog ^ 40 km (see Fig- 
ure 2). Assuming that 5cog = &t (but, see discussion in 
Section 2.1). We find that T > 10'"^ K for absorbers with 
line width 6lw ^ 50 km s~^, not the 6lw > 40 km s~^ 
used by most previous works. This difference is central to 
our analysis. So as a guideline, we established the follow- 
ing categories and criteria for BLAs, using our analysis 
combined with that in the literature: 
A. Probable BLA: There is no definitive test for BLAs, 
but we classify as "probable" any system with well- 
defined 5coG > 40 km s~^ or 6lw > 60 km s~^ confirmed 
by independent measurement of L07 and ourselves and 
with no obvious sign of multiple component structure 
(such as an asymmetric profile) in Lya or higher Lyman 
lines or metal ions (if available). Fifteen systems fall 
within this category. Note that we do not a priori use 
the presence or absence of O VI, N V, or other WHIM 
sign posts as a BLA determinant, although several O VI 
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1494 1496 
wavelength (A) 



_ ^ HS0624+6907 

's 




1612 1614 1616 IBID 

wav clcngth (A) 

Fig. 4. — Examples of three purported Lya lines that do not 
appear in our reduction of the STIS/E140M data. Reported red- 
shifts, f)-values, and column densities listed in L07 are printed in 
the lower left corner of each panel, and the resulting line profiles are 
overplotted. Most of these "missing" lines are associated with very 
narrow bad-pixel regions where the flux drops to negative values 
in our version of the reduced data. 



HE0226-4410 STIS/El lOM 




waTClcngth (A) 

Fig. 5. — Potential "echelle ripple" behavior at the long- 
wavelength end of the STIS/E140M detector. Four Ly« absorbers 
are claimed (Lehner et al. 2006) but only the reddest (z = 0.40027) 
is confirmed with a Ly/3 counterpart. While we view the other 
three absorbers with some skepticism, we find little or no simi- 
lar "echelle ripple" behavior in other datasets and retain them as 
possible BLAs in our accounting. 



detections are present in this sample (see below). 

B. Possible BLA: Less likely than category A, the 
"possible" category includes systems with 40 < 5lw < 
60 km s^^. This category also holds intermediate cases, 
where an absorber shows ambiguous component struc- 
ture or a high degree of uncertainty as to line width 
and/or continuum fit, but a broad H I system compo- 
nent cannot be reasonably ruled out. Forty-eight H I 
absorbers fall in this category. 

C. Non-BLA: Nearly half (56/119) of the H I lines with 
5 > 40 km s^i pubhshed by L07 or DS08 are most likely 
not BLAs by our assessment for a variety of reasons: 



measured 5coG < 40 km s^^ or 6lw < 40 km s^^; prob- 
able alternate line identification; obvious narrow com- 
ponent structure; or simply not detected as absorption 
features in our reduction of the data (i.e., we do not con- 
firm the extraction and/or analysis done by or referenced 
in L07, see Figure 4). The last group, encompassing nine 
purported absorbers, occurs primarily in the spectrum of 
HS 0624+6907 (Aracil et al. 2006a,b). In the course of 
our independent analysis, we identified a few BLA can- 
didates that were not listed in L07 or other independent 
literature sources. In most cases, these lines turn out to 
be something other than Lya and thus fall into category 
C, as they are not confirmed by two independent sources. 

Because the 7 individual sightlines were analyzed orig- 
inally by different authors, somewhat differing analysis 
procedures were used and then adopted by L07. Each is 
briefly discussed below, with notes that bear on the BLA 
identification process. 

HE 0226— 4110 — Originally analyzed by Lehner et al. 
(2006), the spectrum is shown in that paper. No COG 
measurements are reported; instead some 6-values are 
obtained by simultaneous fitting of several Lyman-series 
lines (see discussion below). Longward of 1690 A, the 
continuum in this spectrum appears to "ripple" in both 
our extracted spectrum (Figure 5) that of Lehner et al. . 
This gives rise to four reported BLA candidates (one with 
^LW ~ 150 km s~^) in only Az = 0.007. The periodic na- 
ture of these features are suspicious, especially given that 
these features occur the highest spectral orders which cut 
across many detecor rows before extraction. However, 
no similar features appear in the highest spectral order 
(1690 < A < I710A) for AGN where ZLya > ^agn- 

HS 0624+6907— Originally analyzed by Aracil et al. 
(2006a, b), no spectra were published. These authors do 
not mention obtaining COG 5- values except for a Lya ab- 
sorption complex at 2: ~ 0.06. Therefore, we assume that 
all 6-values quoted by L07 for this spectrum are either 
Lya line width values or simultaneous line width fits to 
two or more Lyman lines as described for HE 0226— 4110 
above. Several of the BLAs reported by L07 in this sight 
line are at or near the locations of detector artifacts; we 
see no believable absorption at these locations based on 
our own reduction of the data. The reported lines may 
be a result of smoothing the data over these bad pixels. 
Longward of 1580 A in this spectrum our extracted data 
does not match even the presence of some of the absorp- 
tion lines listed by L07 (Figure 4). These unconfirmed 
lines are listed in Table 2 as non-BLAs for completeness, 
but are not otherwise analyzed. 

PG 1116+215 — Originally analyzed by Sembach et al. 
(2004), the spectrum is shown in that paper. Sembach 
et al. (2004) use COG 6-values where available, although 
none apply to the BLAs in this sight line. 

PG 1259+593 — This sight line was originally analyzed 
by Richter et al. (2004) and the spectrum is shown in that 
paper. Richter et al. list several COG 6-values for BLAs 
and we add several more. 

PKS 0405-123— Originally analyzed by Wilhger et al. 
(2006, W06) and reanalyzed by L07. Neither W06 nor 
L07 report COG measurements for this sightline. The 
spectrum is shown in Williger et al. Despite mentioning 
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half or less of the L07 reported values; two other sys- 
tems (z = 0.17876 and 0.25861) have similar difFerenees 
between reported values but, with large 6cOG error bars 
and so are not inconsistent with the large line widths re- 
ported by L07. Based upon these 10 systems only, we 
find significant differences between COG &- values and si- 
multaneous line- fit &- values in roughly half the cases, all 
in the sense that the COG 6- values are significantly less 
(35-50% the amount). When higher S/N FUV spectra 
from COS become available, tests using these two tech- 
niques should be made to determine which method is the 
most reliable as a function of line width and S/N. 

3. RESULTS 

All BLA candidates drawn from the catalogs of L07 
and DS08 were independently analyzed by both Dan- 
forth and Stocke based on the criteria above. In Table 2, 
we list the Probable (A), Possible (B), and non-BLA 
(C) absorbers based upon our analyses. Sight line name 
and absorber redshift are in the first two columns. In- 
dependent 6-values and column densities are listed from 
L07 (or similar source, columns 3-5) and DS08/this work 
(columns 6-8). COG measurements of 6hi, A^hi are given 
where possible, otherwise the unweighted mean of the 
two independent Lya measurements is taken as consen- 
sus b, Nm values along with the method used ('LW for 
Lya line width; 'COG' for curve of growth) in columns 
9 and 10. Brief notes are given for each absorber, with 
additional details presented in an Appendix for many ab- 
sorbers in column 11. In total, we find 15 absorbers in 
group A (Probable BLAs) and 48 in group B (Possible 
BLAs) and 56 non-BLAs out of 119 candidate systems 
with 6 > 40 km s~^ as reported by L07, DS08, or sim- 
ilar literature source. We show the Probable BLAs in 
Figure 6 and a selection of Possible BLAs in Figure 7. 

TABLE 2 

Probable, Possible, and non-BLA Absorbers 



Sight Line Zabs h log A^i Src'' b2 log N2 Src'^ 60 log iVo Absorber J 

(km s-l) (cm-2) 1 (km s'^) (cm-^) 2 (km s"!) (cm-^) 



Probable BLA Detections 



HE 0226-4110 


0.40274 


46 + 4 


14.13 + 0.04 


L07 


49 COG 


14.16 + 0.11 


DS08 


49 COG 


14.15 


hya,l3 detections 


PKS 0405-123 


0.08139 


54 + 4 


13.79 + 0.02 


L07 


53 + 2 


13.76 + 0.02 


DS08 


53 LW 


13.77 


sec appendix 


HS 06244-6907 


0.09023 


76 + 14 


13.29 + 0.08 


L07 


100 ; 


13.34 + 0.04 


this 


90: LW 


13.31 


very weak and broac 


HS 0624-1-6907 


0.33976 


42 ±1 


14.45 + 0.03 


L07 


41 COG 


14 46+010 


DS08 


41 COG 


14.45 


OVI, see appendix 


PG 1116-1-215 


0.04125 


105 ± 18 


13 25+° " 


L07 


72 + 6 


10+0. 06 

10.10_Q gg 


DS08 


90: LW 


13.19 


very weak 


PG 1116-1-215 


0.06244 


77 + 9 


1, -,0+0.07 

13.i»_o.o6 


L07 


65 + 7 


13.10 + 0.06 


this 


71 LW 


13.14 


very weak 


PG 1116+215 


0.09279 


133 ± 17 


-,3 og + 0.09 


L07 


83 : 


-,0 n,+0.04 

10.01_Q ,j3 


this 


100: LW 


13.35 


sec appendix 


PG 1116+215 


0.13370 


84 + 10 


IT 97+0.08 


L07 


79 + 6 


13.20 + 0.04 


DS08 


82 LW 


13.24 


OVI, sec appendix 


PG 1259+593 


0.00229 


44^^ COG 


13.61 + 0.06 


R04 


61 + 5 


13.83 + 0.24 


this 


44 COG 


13.61 


sec appendix 


PG 1259+593 


0.14852 


42 + 2 


13.91 + 0.06 


L07 


57 : COG 


13.85 + 0.01 


this 


42 LW* 


13.88 


sec appendix 


PG 1259+593 


0.15136 


65 + 6 


13.32 + 0.09 


L07 


59 + 5 


13.21 + 0.05 


DS08 


62 LW 


13.26 




H 1821+643 


0.11133 


88 + 14 


12.951°;!^ 


L07 


51 : 


12.97 + 0.03 


this 


51: LW 


12.96 


see appendix 


H 1821+643 


0.12147 


40ttt 


14.04 + 0.36 


L07 


76 + 5 


13.79 + 0.06 


this 


70: LW 


13.80 


OVI, see appendix 


H 1821+643 


0.21326 


43 + 2 


14 41+0.04 
-'-^•^-'--0.04 


L07 


42l^ COG 


i'l.+U_Q Qg 


DS08 


42 COG 


14.40 


OVI, sec appendix 


H 1821+643 


0.26658 


45 + 2 


10 54+0.03 


L07 


45 + 1 


13.52 + 0.01 


DS08 


45 LW 


13.58 


OVI, see appendix 


Possible BLA Detections 


HE 0226-4110 


0.06083 


45 + 1 


14.65 ± 0.02 


L07 


46 + 15 


14.86 + 0.15 


this 


45 LW 


14.75 


see appendix 


HE 0226-4110 


0.09220 


40 + 18 


12.94 + 0.11 


L07 


47 + 9 


1, i,+0.08 

iCS.il g gg 


this 


44 LW 


13.03 


weak, profile uncertc 


HE 0226-4110 


0.15175 


49 + 7 


13.42 + 0.05 


L07 


51 + 5 


13.39 + 0.03 


this 


50 LW 


13.41 


asymmetric 


HE 0226-4110 


0.16339 


46 + 2 


14.36 + 0.04 


L07 


39 + 7 


14.34 + 0.03 


this 


42 LW 


14.35 


see appendix 



that the reanalysis was required due to the large number 
of suspect BLA candidates, L07 5-values from simulta- 
neous line fits are similar to the Lya-only line- widths re- 
ported by Williger et al. (2006). We report several COG 
6-values for this sight line and, similar to many other 
absorption systems we generally find 6coG < 6lw (see 
below). 

H I82I+643 For spectral analysis, L07 refer to Sem- 
bach et al. In prep, which does not appear to have been 
published yet. There is no mention of COG measure- 
ments in the L07 description, so we have assumed that 
the reported values are 6lw- 

PG 0953+415 — L07 report these BLAs based upon an 

unpublished Tripp et al. analysis. We have assumed that 
all quoted 6-values are based on Lya line width measure- 
ments since our own 6lw values are similar. 

One of the procedural differences between L07 and the 
current work is that where more than one Lyman line 
is detected in absorption in these spectra, we have used 
a 6-value determined by the COG method whereas L07 
uses a simultaneous line fit of all available lines. In gen- 
eral, the 6coG values are significantly smaller than those 
obtained by the simultaneous line fitting method. To 
evaluate this difference, we have used the 10 absorbers 
in the PKS 0405— 123 sight line analyzed independently 
by ourselves. L07, and W06. For all ten systems the 
W06 (6lw) measurements agree to within the errors with 
the L07 (simultaneous line-fits) 6-values. Half of these 
measurements also show excellent agreement with 6coG • 
However, the other five systems show 6coG values sig- 
nificantly less than values obtained by the simultaneous 
line fit method (for details see Table 2) and these dis- 
agreements are for the largest 6-values reported by L07 
in this sightline. Three systems {z = 0.09659, 0.16121, 
and 0.32500) have COG 6-values with tight error bars 
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TABLE 2 

Probable, Possible, and non-BLA Absorbers 



HE 0226-4110 


0. 


,18619 


54 + 16 


13.26 ± 0.08 


L07 


37 ±8 


13.02 ±0.07 


this 


45: LW 


13, 


.14 


possibly two compor 


HE 0226-4110 


0. 


,20700 


97 : 


13.31 ± 0.39 


805 


blend 


< 15.20 


this 


97: LW 


13, 


.31 


OVI, NeVHI, see ap 


HE 0226-4110 


0. 


,30930 


44 ±2 


14.26 ± 0.03 


L07 


COG 


U.28+°i° 


this 


40:COG 


14, 


.27 


see appendix 


HE 0226-4110 


0. 


,39641 


63 ±23 


13.59 ±0.10 


L07 


57± 10 


13.54 ±0.10 


this 


60 LW 


13, 


.57 


see appendix 


HE 0226— 4110 


0. 


,39890 


152 : 


13.50 ± 0.16 


L07 


88 ± 20 


13.46 ± 0.10 


this 


100: LW 


13, 


.50 


see appendix 


WTT noofi y1 1 1 n 
jlCi UxzO — 4tllU 


U. 




Dl It ZD 


1 Q QQ + n 1 1 


T n7 


1 4 It ID 


1 Q OQ + n 1 1 

io.zy It u.il 


tnis 


7rv 


1 Q 




see appendix 


PKS n4n'i— 123 


0. 


,03196 


54 ± 16 


1 Q Qo + n n« 

XO.tJO _1_ U.L»0 


L07 


57 ± 8 


1 Q QQ + n nfi 

XO.tJO _1_ \J.\J\J 


DS08 


56 LW 


13 


.33 




PKS 0405-123 


0. 


,07523 


48 ±20 


13.05 ±0.11 


L07 


37± 12 


12.81 ±0.13 


this 


43 LW 


12, 


.93 


see appendix 


PKS 0405—123 


0. 


09659 


70 ± 20 


13 90 + 18 


L07 


36 + 4 COG 


14 64+° - 12 

-0.11 


DS08 


70: LW 


13, 


.90 


v_y V x^ ov^V' Op^^^xxvii-A. 




n 
u. 


13102 


□z zn o 




T 07 


oo It u 


1 Q + n f\A 
±0.00 It U.U4 


f liic 

tnib 


t>7 TW 


1 'i 


41 


see Appendix 


PKS n4f)5— 123 


0. 


,13377 


43 ± 8 


13 34 + (16 


L07 


45 + 7 


13 22 + 07 


this 


44 LW 


13, 


.28 


OV^V^ .IL L/L^^XXVAl-A. 


PKS 0405-123 


0. 


,16678 


75 ±7 


13.91 ± 0.04 


L07 


77 ±7 


13.84 ± 0.05 


this 


76 LW 


13, 


.88 


poss. OVI, see appei 


PKS 0405-123 


0. 


,17876 


55 ±7 


13.61 ± 0.04 


L07 


18t^* COG 


13.92i;;;2? 


DS08 


55 LW 


13, 


.62 


asymetric, see apper 


PKS 0405-123 


0. 


,18269 


48 ± 2 COG 


15.07 ± 0.09 


W06 


49+10 COG 


14 90+',l if; 


DS08 


49 COG 


14 


.86 


OVI, see appendix 


PKS 0405—123 


0. 


,19086 


44 ± 16 


13.17 ± 0.09 


L07 


38 ± 9 


1 o ni +0.10 
13.01Io.o6 


this 


41 LW 


13 


,09 


sec Appendix 


PKS 0405-123 


0, 


,24513 


54 ± 24 


13.23 ± 0.11 


L07 


30 ± 5 


13.01 ± 0.09 


this 


42: LW 


13, 


,12 


poor agreement, see 


1>KS 0405-123 


0, 


,25861 


40 ±9 


13.37 ±0.07 


L07 


39 ±4 


1 3 qfi+O-ls 


this 


40 LW 


13, 


,36 


bcoG = 2ll?g; woe 


PKS 0405-123 


0, 


,29523 


47± 13 


13.33 ±0.08 


L07 


39 + 8 


13.07 ±0.08 


this 


43 LW 


13, 


.20 


PKS 0405-123 


0, 


,29904 


49 ± 23 


13.26 ±0.12 


L07 


54 + 10 


13.14 ±0.15 


this 


52: LW 


13, 


.20 


noisy, components? 


PKS 0405—123 


0, 


,35092 


38 ± 2 


14.25 ± 0.03 


L07 


56 ± 15 COG 


14.04 ± 0.08 


DS08 


40: LW 


14, 


.25 


sec appendix 


PKb 0405—123 


0. 


40886 


40 ± 2 


14.38 ± 0.03 


L07 


40_!_4 COG 


14.32 ± 0.05 


DS08 


40:COG 


14, 


,35 


sec appendix 


Hb 0024+0907 


0. 


05437 


60 ± 19 


13.09 ± O.li 


L07 


45 ± 5 


i^.»»_(, jj 


this 


52: LW 


12 


,99 


sec appendix 


HS 0fi24+fiQn7 


0. 


05515 


84 ± 31 


1 3 fiS + n 1 7 


L07 


87 ± 13 


1 q «q + f) 06 


this 


85- T,W 

. J_J V V 


13 


,66 




HS 0624+6907 


0, 


13597 


57 ± 11 


13.33 ± 0.10 


L07 


52 ± 12 


13.16 ± 0.10 


this 


55 LW 


13, 


,25 


OVI, see appendix 


HS 0624+6907 


0, 


,21323 


45 ± 6 


13.22 + 0.05 


L07 


40 ± 7 


13.04 ±0.08 


DS08 


43 LW 


13, 


,13 


poss. Siiii, cm 


HS 0624+6907 


0, 


,26856 


51 ± 7 


13.03 ± 0.05 


L07 


46 ±8 


13.00 ±0.07 


this 


49 LW 


13, 


.02 


weak 


HS 0624+6907 


0, 


,29661 


52 ±3 


13.54 + 0.02 


L07 


45 


13.31 + 0.04 


DS08 


48 LW 


13, 


,43 




HS 0624+6907 


0, 


,30994 


66 ± 12 


13.61 ± 0.10 


L07 


38 ± 9 


12.90 ± 0.11 


DS08 


52: LW 


13, 


.25 


poss. OVI, see appei 


HS 0624+6907 


0, 


,31790 


34 ±4 


13.37 ±0.04 


L07 


51 ±7 


13.381°;°^ 


DS08 


43 LW 


13, 


.38 


OVI, see appendix 


PG 0953+415 


0. 


,05879 




13.41 ±0.16 


L07 


not measured 


not measured 




63: LW 


13, 


.37 


see appendix 


PG 0953+415 


0. 


,17985 




13.27 ±0.07 


L07 


48 ±2 


13 20+°"^ 


DS08 


48 LW 


13 


.24 


see appendix 


PG 0953+415 


0, 


,19126 




13.08 ±0.55 


L07 


66 : 


13.50 ±0.04 


this 


40: LW 


13, 


,30 


sec appendix 


PG 0953+415 


0, 


,20104 




13.16 ±0.16 


L07 


44 : 


12.63 ±0.19 


this 


57: LW 


12, 


.90 


weak, uncertain coni 


PG 1116+215 


0, 


,01638 


49 ± 5 


13.39 ±0.06 


L07 


54 ±4 


13.44 ±0.03 


DS08 


51 LW 


13, 


,41 


asymmetric 


PG 1116+215 


0, 


,06072 


55 + 6 


13.281" (Js 


L07 


53 


13.21 ±0.04 


DS08 


54 LW 


13, 


,24 




PG 1116+215 


0. 


,08587 


52 ± 14 


12.90to.l3 


L07 


55 ± 11 


12.88 ± 0.08 


this 


53 LW 


12 


.89 


see appendix 


PG 1259+593 


0. 


,19573 


43 ± 7 


13.07 ±0.14 


R04 


46± 19 


13.05 ±0.15 


DS08 


44 LW 


13, 


.06 


see appendix 


PG 1259+593 


0. 


,41786 


51 ±4 


13.25 ±0.08 


L07 


52 ±24 


13.23 ±0.18 


DS08 


51: LW 


13, 


.24 


asymmetric 


H 1821+643 


0. 


,02642 


49 ±6 




L07 


44 ±3 


13.19 ±0.03 


DS08 


47 LW 


13, 


.23 


asymmetric 


H 1821+643 


0. 


,14754 


45 ±3 


lo CI +0.03 
10.01_o Q3 


L07 


42 ±2 


13.43 ± 0.02 


this 


44 LW 


13 


.47 


see appendix 


H 1821+643 


0, 


16352 


52 ± 6 


iq 1 7+0.06 
lO.ll _o Q7 


L07 


56 ±4 


13.11 ±0.04 


DS08 


54 LW 


13, 


,14 


weaJf, slightly asymr 


H 1821+643 


0, 


,18047 


51 ± 7 


lCi.i4_g Qg 


L07 


52 ±8 


13.05 ±0.07 


DS08 


51 LW 


13, 


.10 




H 1821+643 


0. 


,22616 


55 ±4 


i3.5iini 


L07 


56 ±2 


13.41 ± 0.07 


this 


56 LW 


13, 


.46 


offset OVI, see appe: 


H 1821+643 


0. 


,25814 


60 ±9 


13.38l°;l°3 


L07 


56 : 


13.26 ± 0.04 


this 


58: LW 


13, 


.32 


see appendix 



Not BLAs 



HE 0226-4110 


0, 


,01216 


not detected 


not detected 




41 + 7 


1 Q nc+0.08 


DS08 








poor S/N, not confir 


HE 0226-4110 


0, 


,02679 


42 ± 11 


13.22 ±0.08 


L07 


30 ± 4 


13 17+" 12 


DS08 


36 LW 


13, 


.19 


Lehner et al. 2006:6] 


HE 0226-4110 


0, 


,23009 


68 ±8 


13.69 ± 0.04 


L07 


not detected 


not detected 










no absorption seen 


HE 0226-4110 


0. 


,22102 


34± 18 


12.99 ±0.12 


L07 


44 ±9 


13.06in^ 


this 


39 LW 


13, 


.02 


weaJf 


HE 0226-4110 


0. 


,38420 


62 ±7 


13.91 ± 0.04 


L07 


31 : COG 


13.08lg;li 


this 


31: COG 


13, 


.49 


components, see app 


PKS 0405-123 


0. 


,05896 


77 ±27 


13.34 ±0.10 


L07 


20,20 


12.97 ±0.08 


this 


20 LW 


12, 


.97 


see appendix 


PKS 0405-123 


0. 


,07218 


45 ±14 


13.09 ± 0.08 


L07 


45 : 


13.01 ±0.06 


this 


45 LW 


13, 


.05 


see appendix 


PKS 0405-123 


0. 


,10298 


87± 19 


13.40 ± 0.07 


L07 


54 : 


1 o oc: + 0,04 

10.00_Q Q3 


this 


70: LW 


13, 


.38 


see appendix 


PKS 0405-123 


0. 


,10419 


not detected 


not detected 




48 ±5 


12.98 ± 0.08 


DS08 








see appendix 


PKS 0405-123 


0. 


,13646 


54± 11 


13.34 ± 0.06 


L07 


11,19 


12.96 ±0.07 


this 


19: LW 


12 


.96 


see appendix 


PKS 0405-123 


0. 


,13924 


35 ± 11 


13.02 + 0.11 


W06 


53 ± 2 


,o ,0+0.06 
-'-■'•-'-■'-0.05 


DS08 


41 LW 


13 


,08 


sec appendix 


PKS 0405-123 


0, 


,15304 


46 + 3 


13.80 + 0.03 


L07 


50 : 


13.76 ±0.03 


this 


48: LW 


13, 


.78 


sec appendix 


PKS 0405-123 


0, 


,16121 


54 + 8 


13.71 + 0.04 


L07 


27± 14 COG 


13.62lf;;?i 


DS08 


27 COG 


13 


.66 


see appendix 


PKS 0405-123 


0, 


16714 


30 ± 1 


16.27 ±0.13 


L07 


40tJ.2 COG 


.„+0.40 

io.+ r_(, 4g 


DS08 








OVI, components, sc 


PKS 0405-123 


0, 


,24057 


57± 19 


13.27 ±0.09 


L07 


58 ±6 


-|o 97+0.08 

iO.Z(_Q Q7 


this 


58 LW 


13, 


,27 


marginal feature, ar1 


PKS 0405-123 


0, 


,28838 


52± 19 


13.32 ±0.10 


L07 


30 : 


13.18 ±0.10 


this 


30: LW 


13, 


.10 


components; W06:6i 


PKS 0405-123 


0. 


,32500 


66± 13 


13.55 ± 0.06 


L07 


60 : 


lCS.40_o Qg 


this 


63: LW 


13, 


.50 


see appendix 


PKS 0405-123 


0. 


,34234 


42 ±13 


13.39 ± 0.08 


L07 


35 ±4 


,Q 97+0.07 
lcS.Z/_(, 05 


this 


38 LW 


13, 


.33 


asymetric 


PKS 0405-123 


0. 


,36150 


44± 10 


13.71 ±0.10 


L07 


25 ± 3 COG 


15.18 ±0.13 


this 


30: COG 


15, 


.00 


see appendix 


HS 0624+6907 


0. 


,04116 


41 ±3 


13.33 ± 0.03 


L07 


not measured 


not measured 










SIIIz = 0.0635, see a 


HS 0624+6907 


0. 


,05483 


35 : 


14.50 : 


L07 


45lf COG 


14.281^1^ 


DS08 


35: LW 


14, 


.28 


see appendix 


HS 0624+6907 


0. 


,06346 


48 ±8 


14.46 ± 0.30 


L07 


33t? COG 


15 25+°-l° 


DS08 


33 COG 


15, 


.25 


OVI, see appendix 
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TABLE 2 

Probable, Possible, and non-BLA Absorbers 



HS 0624+6907 


0.19979 


17 ±2 


13.24 ± 0.05 


L07 


55 + 10 


13.58 + 0.03 


DS08 


17 LW 


13.41 


poss. OVI; compone 


HS 0624+6907 


0.21990 


60 ±9 


13.39 ± 0.05 


L07 


not detected 


not detected 








not detected 


HS 0624+6907 


0.23231 


44 + 8 


13.33 ± 0.08 


L07 


not detected 


not detected 








not detected 


HS 0624+6907 


0.28017 


43 + 2 


14.32 ± 0.02 


L07 


35^5 COG 




DS08 


35 COG 


14.37 




HS 0624+6907 


0.29531 


42 + 2 


13.80 ± 0.02 


L07 


34 


13.62 + 0.03 


DS08 


38 LW 


13.71 




HS 0624+6907 


0.31045 


62 ± 40 


13.43 ± 0.33 


L07 


not detected 


not detected 








see appendix 


rib 0624+6907 


0.31088 


51 ± 28 


IS.IS ± OAS 


L07 


not detected 


not detected 








see appendix 




n '^i 9sn 

U.O ±ZrOU 


54 + 9 




L07 












iioi ucteeicti 


HS 0624+6907 


0.31326 


55 + 11 


13.62 + 0.10 


L07 


not detected 


not detected 








not detected 


HS 0624+6907 


0.32089 


31 + 1 


13.97 ± 0.02 


L07 


44j:J| COG 


13.80+^?? 

—0.11 


DS08 


31 LW 


13.89 


see appendix 


HS 0624+6907 


0.32724 


69 ± 16 


13.73 ± 0.32 


L07 


not detected 


not detected 








not detected 


HS 0624+6907 


0.32772 


115 + 62 


13.61 ± 0.43 


L07 


not detected 


not detected 








not detected 


PG 0953+415 


0.02336 




13.21 ± 0.08 


L07 


20 : 


12 68+°-l? 


this 


38: LW 


12.95 


discrepant measuren 


F(j 0953+415 


0.04382 


not detectea 


not detected 




47 + 8 


1 Q qQ+0.03 
io.oy_g 04 


DS08 


47 LW 


13.06 


UVl2 = 0.22974; see 


PG 0953+415 


0.12784 


44 : 


12.83 ± 0.35 


L07 


not detected 


not detected 


this 






not detected, see ap] 


PG 0953+415 


0.19361 


40 + 2 


13.94 + 0.02 


L07 


37 + 2 


13.82 + 0.03 


this 


39 COG 


14.15 


see appendix 


PG 0953+415 


0.20006 




13.24 + 0.09 


L07 


32 ± 7 


12.91 + 0.11 


DS08 


49: LW 


13.07 


continuum uncertain 


PG 1116+215 


0.02841 


31 ± 1 


13.80 + 0.02 


L07 


43 : 


1 3 71 +0-21 
' -^-0.09 


DS08 


35 LW 


13.78 




PG 1116+215 


0.05904 


21, 30 


??? 


804 


45 : 


13.53 + 0.02 


DS08 


21,30 LW 


13.53 


components 


Pvj Izoy+oyo 


U.U4oUO 


48 + 12 


lo.oo ± U.zl 


L07 


Qn + -'-0 r^r\r^ 

3U_^ UULt 


-ic: C--I+0.28 
i0.01_Q 25 


JJbOo 




15.55 


OVI, components, sc 




u.uDyoi 


not detected 


not detected 




0^ JZ ^ 


13.33 + 0.04 








not detected by P04 


PG 1259+593 


0.08041 


42 + 4 


12.97 + 0.10 


L07 


26 + 2 


12.86 


this 


34: LW 


12.92 


not detected by R,04 


PG 1259+593 


0.17891 


99 ± 9 


13.29 + 0.10 


L07 


not measured 


not measured 








not detected 


PG 1259+593 


0.21136 


not measured 


not measured 




47 + 3 


13.39 + 0.04 


DS08 


47 LW 


13.39 


see appendix 


PG 1259+593 


0.22861 


40 + 3 


13.47 + 0.05 


L07 


34 + 2 


13.42 + 0.03 


DS08 


37 LW 


13.14 




PG 1259+593 


0.24126 


89 + 7 


13.41 + 0.09 


L07 


23 + 4 


12.86 + 0.08 


this 


55 LW 


13.14 


ambiguous compone 


PG 1259+593 


0.25971 


41 ± 5 


13.84 + 0.12 


L07 


29+^ COG 


13.98 + 0.07 


DS08 


29 COG 


13.91 


OVI, obvious compo 


PG 1259+593 


0.30434 


65 + 1 


13.76 + 0.14 


L07 


20tl OOG 


-iq O1+0.07 


DS08 


20 COG 


13.81 


obvious components 


PG 1259+593 


0.31978 


74 + 8 


13.98 ± 0.06 


L07 


31 + 3 COG 


14.07lH^ 


DS08 


31 COG 


14.07 


OVI; listed as blend 


PG 1259+593 


0.32478 


46 + 10 


13.24 + 0.15 


L07 


20 


12.93l°;lt 


DS08 


33: LW 


13.10 


OVI; continuum unc 


PG 1259+593 


0.43569 


44 + 4 


14.22 + 0.10 


R04 










see appendix 


H 1821+643 


0.12221 


42 + 4 


13.19 + 0.06 


L07 


34 + 5 


iq iq + 0.03 

10.10_g Q4 


DS08 


38 LW 


13.16 




H 1821+643 


0.19176 


not detected 


not detected 




47 + 7 


12.72 + 0.07 


DS08 






marginal; not confiri 


H 1821+643 


0.22480 


not detected 


not detected 




40+1 OOG 


15 41+0-25 


DS08 


40 COG 


15.41 


OVI; components, se 



Measurement sources: L07 (Lehner et aL 2007 and sources 
therein); DS08 (Danforth & Shull 2008); W06 (Williger et al. 
2006); 805 (Savage et al. 2005); S04 (Sembach et al. 2004); R04 
(Richter et al. 2004); this (this work). 

3.1. Broad Lya Absorber Frequency 

In measuring metal-ion absorption linos, DS08 employ 
a 4(7 detection limit, with equivalent width 
^ jSL) jX/R) 

- {S/N)x ' 
where the instrumental resolving power R = A/ A A and 
SL is the significance level in standard deviations. How- 
ever, both DS08 and other studies locate Lya lines inter- 
actively, which does not follow the strict Wmin argument 
used for lines of metal ions. Determining the significance 
level of the Lya lines reported by DS08 based a posteri- 
ori on the observed equivalent widths and data S/N, we 
find SL > 10—15 for the weakest detections. Weaker 
features arc certainly visible in the data, but they can 
reasonably be explained as fixed-pattern noise or other 
instrumental features. Indeed, DS08 use SL > 10 when 
determining the rcdshift pathlcngth AzLya in their ab- 
sorber frequencies and c;osnioIogical calculations. 

Tripp ct al. (2008) point out that equation 5 is strictly 
valid only for unresolved features. Since BLAs are several 
times wider than the instrumental resolution element for 
both STIS and FUSE data, equation 5 is even less accu- 
rate. However, we argue that, even if the SL is not rigor- 
ously correct, it is still relatively correct, and it gives us 



a basis for equal comparison between lines. Broader lines 
will be shallower for the same equivalent width and thus 
each pixel will show less contrast from the continuum 
at full spectral resolution. Smoothing over the number 
of pixels equal to the line width (something the human 
brain does naturally) to first order yields no change in 
sensitivity for lines of the same equivalent width but dif- 
ferent fe- values. 

Previous BLA studies used the quantity log(A^/6) as 
a detection criterion, reasoning that broader lines re- 
quire a higher column density (and hence equivalent 
width) to reach the same line center optical depth. In 
particular, Richter et al. (2006) used log{N/b) > 11.3 
[for N in cm~^ and b in km s~^] as their detection 
threshold. For a constant column density, log {N/b) 
changes rapidly for narrow lines, but much more slowly 
at 6 > 40 km s~^. For example, log{N/b) > 11.3 corre- 
sponds to log Ni^ya > 12.9 at 6 = 40 km s~^ but only 
0.2 dex (60%) higher at 6 = 65 km s~^. We argue that 
setting a detection threshold based purely on equivalent 
width, Wlyq {S/N)/AX > 4, is roughly equivalent to 
setting one in (N/b) for BLAs. It should be noted that 
^ 75% of the Possible and Probable BLAs in our sample 
show log(iV/6) > 11.3, and all have log(iV/6) > 11.0. 

Total absorption pathlength Az depends on the red- 
shift of the background AGN and the wavelength cov- 
erage of the UV spectrograph. We follow a procedure 
identical to that discussed in DS08; the local S/N in the 
data is defined as the mean flux divided by the standard 
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b=49 COG 
logN=14.15" 
= 0.40274 ^LC^'" probable BLA 

















CO 
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d 














b = 53 LWh 


PKS0405-123 




logN=13.77 : 


z = 0. 08139 




probable BLA = 
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05 

o 
o 




HS0624+6907 
z = 0. 09023 


b=90 : LW- 
logN=13.31 ; 
probable BLA : 



1704.0 1704.5 1705.0 1705.5 1706.0 
Wavelength 



1313.0 1313.5 1314.0 1314.5 1315.0 
Wavelength (K) 



1315.5 1316.0 



1324.0 1324.5 1325.0 1325.5 1326.0 1326.5 
Wavelength (A) 




b = 41 COG 
logN=14.45" 
probable BLA 



1627.5 162B.0 162B.5 1629.0 1629.5 1630.0 
Wavelength (A) 
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"PG1116+215 
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z = 0. 04125 


probable BLA 
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1265.5 1266.0 1266.5 1267.0 




Wavelength (A) 



PG1116+215 
z = 0. 06244 



b = 71 LW: 
logN=13.14 
probable BLA 



1290.5 1291.0 1291.5 1292.0 1292.5 1293.0 
Wavelength (A) 



PG1116+215 
z = 0. 09279 



1327.0 1327.5 1328.0 1328.5 1329.0 
Wavelength (A) 
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logN=13.35 
probable BLA - 

329.5 1330.0 




PG1116+215 
z = 0. 13370 

1377.0 1377.5 



1378.0 1378.5 
Wavelength (A) 



b=82 LW 
logN=13.24 
probable BLA 

379.0 1379.5 



2":PG1259 + 593 
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1.5 - 



0.5 - 




b=44 COG: 
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probable BLA 

1218.0 1218.5 1219.0 1219.5 1220.0 
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b=42 LW* : 
logN = 13.88T 
probable BLA : 
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= 0.15136 
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logN=13.26 
probable BLA 



_H1821+643 
z = 0. 11133 



b=51 : LW 
logN=12.96 
probable BLA 



1395.5 1396.0 1396.5 
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Wavelength (A) 
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CM 
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d 
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d 
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.1... 


d 






\ 

\ 






:H1821+643 
: z = 0. 12147 
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logN=13.80 ■ 
probable BLA 


1362.0 1362.5 


1363.0 
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z = 0. 21326 
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logN=14.40 
probable BLA 



Wavelength (A) 



1473.5 1474.0 1474.5 1475.0 1475.5 1476.0 1476.5 
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tQ 
Cvl 
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/ b=45 LW 


"H1821+643 




logN=13.58 ■ 


z = 0. 26658 




probable BLA 


1538.5 1539.0 


1539.5 
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Wavelength (A) 



Fig. 6. — Spectra of all 15 Probable BLA systems illustrate the diversity of line profiles. The sight line and absorber redshift are listed on 
the left side of each panel while the consensus 6 (km s~^) and log N values from Table 2 are listed on the right. Dotted line shows continuum 
fit to the data, and dashed curve shows consensus absorption profile. Adjacent line detections are indicated with vertical ticks and labeled 



with ion and redshift. Lines identified without redshift denote Galactic interstellar absorption. Flux is in units of 10 
See Table 2 and Appendix for more details of individual systems. 



erg cm 



deviation in continuum regions after the data have been 
smoothed to the resolution element. The S/N{X) is then 
modified by setting regions with strong IGM, Galactic, 
and intrinsic absorption systems and instrumental fea- 
tures equal to zero. Wmia^hyaiz) is calculated from the 
S/N{X) vector. The path'length Az(Wmin) is then the 



sum of pixels where W > Wmin- As in DS08, we omit 
regions within 500 km of the Galaxy and within 1500 
km s~^ of the AGN to eliminate (most) absorbers intrin- 
sic to either the AGN or the Local Group. Cosmolog- 
ically corrected pathlength AX is calculated in an en- 
tirely analogous manner, using dX = (1 + z)^ [i^jn(l + 
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z = 0. 15175 
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logN=13.41 ; 
possible BLA ; 
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b = 55 LW^ 
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logN=13.24 
possible BLA 
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logN=12.89 
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Wavelength (A) 



1436.0 1436.5 



1490.0 1490.5 1491.0 
Wavelength (A) 



1491.5 1492.0 



Fig. 7. — Same as Fig. 6 but for 15 of the 48 Possible BLA systems. Additional Possible BLAs are shown in Figure 10. 



z)'^+ilA] dz. Throughout this paper we assume a flat 
{rim, A) cosmology with Hq = (70 km s~^ Mpc~-^) hyo, 
n,n = 0.261, r^A = 0.716, and fib = 0.0455 /i^ (Spergel 
et al. 2007). The maximum path lengths available in 
each sight line for absorbers of any strength are listed in 
Table 1. The total pathlength surveyed in all seven sight 
lines is Aztot = 2.193 (AXft = 1.773). 
For the numerator of (cW /(iz)BLA, we have several 



options to choose from, depending how much faith we 
place in our BLA designations. The most skeptical view 
accepts only our Probable sample, which with one-sided 
Poisson statistics gives A/b LA = 15j^| and dAf /dz — 7±2. 
A more inclusive view includes both the Probable and 
Possible groups: A^bla = 63l^ and dM /dz = 29 ± 4. 
Since ^ 50% of the BLAs survive the statistical correc- 
tion process described in Section 3.2, we adopt instead an 
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intermediate census: all of the Probable BLAs and 50% 
of the Possible sample: Mbla = 39^q. Given the micer- 
tainties surrounding BLA identification, we believe pure 
Poisson uncertainties are far too optimistic. We adopt 
the skeptical and inclusive censuses as our lower and up- 
per bounds: A^bla = 39 ±24, {dN'/dz)BLA = 18±11 or, 
in comoving coordinates, (dJ\f/dX)BLA = 22 ± 14. De- 
spite the corrections above, uncertainties in pathlength 
are small (< 10%), so we ignore errors in the denomina- 
tor. 

3.2. Overlap between BLAs and Metal-Line Absorbers 

Because both BLAs and highly-ionized metal lines 
(O VI, N V, Ne VIII, etc.) are thought to trace WHIM 
gas, it is instructive to look at the overlap in these two 
samples. O VI has by far the best detection statistics of 
any FUV intergalactic metal line, with ^ 100 detections 
in the low-z IGM surveys (Danforth & ShuU 2008: Tripp 
et al. 2008; Thom & Chen 2008). In Stocke et al. (2007), 
we used log iVovi > 13-2 and log A^ovi < 13.2 as the 
threshhold between a good O VI detection and a reliable 
non-detection based on the distributions of observed de- 
tections and 4(7 upper limits. Using the same threshold, 
we find that four probable BLAs show O VI absorption, 
while six show O VI non-detections; a detection rate of 
~ 40%. For the larger sample of probable-plus-possible 
BLAs, the numbers rise to 8 detections and 33 non- 
detections for an O VI detection rate of ^ 20%. Using 
the same criteria, the large DS08 survey (~ 650 H I sys- 
tems) features 69 O VI detections and 293 non-detections 
(19%). However, this sample mixes broad with narrow 
H I lines. If we instead define a control sample of non- 
BLAs as all DS08 H I systems with 6hi < 40 km s"^ 
(516 systems), there are 47 O VI detections and 245 non- 
detections (16% detection rate, half that of our Probable 
BLA sample). The O VI detection rate is slightly lower 
(14%) if we use the DS08 6lw measurements to define 
the non-BLA control group. 

The coincidence of probable BLAs and O VI detections 
(40%) is several times higher than in the larger, narrow 
Lya absorber sample (~ 15%). This suggests that BLAs 
and O VI are tracing the same material, but the small 
number of systems in the BLA sample reduces the signif- 
icance. It is worth noting that the 40% O VI detection 
rate is a bit higher than the ~ 25% fraction of H I sys- 
tems that show metal absorption in any ion reported by 
DS08. This suggests that BLAs are accurately tracing 
WHIM irrespective of metal enrichment. Unfortunately, 
the detection statistics for other ions are too poor to draw 
any conclusions. 

3.3. The Galaxy-BLA Connection 

One of the strongest potential advantages of detecting 
the WHIM using BLAs is that these absorbers are not 
affected by the metallicity of the gas, so that even metal- 
free WHIM can be detected. We would expect these low- 
metallicity and metal-free absorbers to be found in re- 
gions far from galaxies, unlike the O VI absorbers which 
are typically found within ^ 800 kpc of the nearest L* 
galaxy (Stocke et al. 2006; Wakker & Savage 2009) and 
even closer to sub-L* galaxies (Stocke et al. 2006). Un- 
fortunately, only eight of the BLAs reported here are 
found in sky regions surveyed for galaxy redshifts com- 
plete to L* or below; these are listed in Table 3 in in- 



creasing order of galaxy separation. Dividing this sample 
into O VI detections and non-detections at a consistent 
level of logA^ovi = 13.2 (DS08), we find nearest galaxies 
at 0.75 — 2.9 h^^ Mpc for the O VI non-detections and 
0.06-0.49 ftfo^ Mpc for the dctctions. Therefore, we find 
some evidence that BLAs are tracing WHIM gas more 
remote from galaxies than by using O VI absorption as 
a WHIM tracer. 

TABLE 3 

Galaxy-BLA Relationship in Well-Surveyed Fields 



AGN 




&Hl" 


BLA 


logiVovi'' 


d 






(km s-l) 


class 


(cm-2) 


(Mpc) 


PG 1259-1-593 


0.00229 


44 COG 


A 


13.7:== 


0.06 


PKS 0405-123 


0.16678 


75: LW 


B 


14.1 + 0.2 


0.11 


PKS 0405-123 


0.09659 


70: LW 


B 


13.7 + 0.2 


0.27 


PKS 0405-123 


0.08139 


53 LW 


A 


13.3 + 0.2 


0.49 


PG 1116+215 


0.06072 


54 LW 


B 


< 13.07 


0.75 


PG 1116+215 


0.09279 


100: LW 


A 


< 13.07 


1.3 


PG 1116+215 


0.01635 


51 LW 


B 


< 13.05 


2.0 


PG 1116+215 


0.08587 


53 LW 


B 


< 13.04 


2.9 



^ Consensus b value from Table 2 

^ OVI column density from DS08 
Detection from Riclitcr ct al. (2004) based on nighttime-only 
FUSE data. Formal significance level is low (< 3(t), however 
absorption appears over an exceptionally broad velocity range 
(-110 <v< +220 km s"!). 

There is one slightly controversial absorber that we 
have counted as an O VI detection in the above account- 
ing: the z = 0.00229 Lya absorber toward PG 1259+593 
is 60 kpc from the 14*'*-mag edge-on late-type spiral 
galaxy UGC8146. DS08 report this absorber as an O VI 
non-detection according to their Aa detection threshhold, 
but Richter et al. (2004) report a low-significance, very 
broad O VI detection at A^'ovi = 5 x lO^'' cm^^ based on 
night-only FUSE data. We thus include this as an O VI 
detection. 

4. BLA COSMOLOGY 

Only about half of the baryons can be accounted for 
in the local universe. The Lya forest makes up only 
about 30% of the total predicted baryons at z ~ (Pen- 
ton, Stocke, & ShuU 2004; Lehner et al. 2007; Danforth 
& ShuU 2008) while collapsed structures (stars, galax- 
ies, etc.) make up another ^ 10% (Salucci & Persic 
1999). Much of the remainder is expected to lie in the 
ionized phases of the IGM above 10^ K. The ~ 100 
O VI absorbers observed at z < 0.4 have been used to 
trace the WHIM phase and can account for an additional 
< 10% of the baryons (DS08), but this estimate relies 
on metallicity and ionization-fraction assumptions that 
make the quantity uncertain (Danforth 2009). The O VI 
WHIM surveys require metal enrichment, leaving open 
the possibility that a significant population of metal-poor 
IGM clouds may contribute to the baryon census. The 
strength of BLA surveys is their ability to trace gas at 
T = 10^ — 10^ K independent of chemical enrichment. 
While there is some overlap with O VI WHIM absorbers, 
BLAs open a new window on the cosmic baryon census. 

4.1. Baryon Fraction Traced by Broad HI 

The mass fraction of the local universe traced by broad 
Lya absorbers can be determined by dividing the to- 
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tal hydrogen column density by the total observed path- 
length 

Obla = ( ^^^^] fi^ . (6) 



C Pcrit 



Since the vast majority of IGM hydrogen is ionized, to- 
tal hydrogen column density can be approximated for 
any given absorber by TVh « -^hi//hi(T)- Neutral hy- 
drogen column Nm can be measured directly in most 
cases. However, the hydrogen neutral fraction /hi is de- 
termined by both photoionization from the metagalactic 
radiation field and ionization due to (thermal) electron 
collisions. At 2: « 0, the ionizing background produces 
an H I photoionization rate Th = 3.2+^^ x 10"^'' s~^ as 
derived in ShuU et al. (1999) from populations and ra- 
diative transer calculations of Seyferts, QSOs, and star- 
bursts. For electron impact, the hydrogen ionization rate 
can be approximated as 



Ci{T) = (5.83x10 



„ , 1, exp[-1.58 X loVr] 
""cm'^s"^) ^ ' ' 



1 + 1.58 X lOVT 



(1.3 X 10"^) [ne/10"*^"cm-3j we adopt 



Th = 3.2 X lO"^'^ s-i and an = 4.1 x lO'^^^ cm^ s'^ 



(7) 

where 1.58 x 10^ X = (/ff/fcs) = (13.6 eV)/kB. The 
critical density where coUisional ionization equals pho- 
toionization is then (ne)crit = TH/Ci{T). For borderline 
WHIM temperatures (T = 10^ K), the ionization rate 
is Ci = 3.6 X 10^'' cm"^ s^^ and the critical density is 
('T'e)crit = 8.9 X 10^^ cm^'^. Howcvei, collisional ioniza- 
tion becomes more and more dominant at higher tem- 
peratures, and at log T = 5.5 (6.0), Ci{T) = 1.7(3.1) x 
10"^ cm^ s"^ and (ne)crit = 1.9(1.0) x 10"^ cm'^ corre- 
sponding to overdcnsitics of (5 < 10 at z ^ 0. 

Photoionization is thus an important consideration 
at low temperatures (fcT <C 1 Ryd), but at tempera- 
tures near or above the O VI peak in CIE (log T = 
5.5 - 6.0), we have (ne)crit = few x 10~^ cm~^. At 
those low densities, in photoionization equilibrium, /hi = 
ne {aH)/{TH 

(case-A at 10** K). If He « (ne)crit) collisional ionization 
could double the ionization rate and halve /hi in the for- 
mula above (0.65 x 10~^). A BLA with rig = (?T.e)crit and 
Nhi = 3 X 10^^ cm-2 would have Nh 5 x 10^^ cm'^ 
and line-of-sight dimension Lbla = Nh/uh ^ 5 x 
10^'' cm =1.5 Mpc. An unvirialized cloud this size would 
exhibit a 100 km s~^ broadening of its Lya linewidth due 
to Hubble flow from one side to another. 

We adopt fm{T) values derived from a set of CLOUDY 
simulations (solid curve in Fig. 1) featuring both colli- 
sional and photoionization (log U = —2, typical of the 
low-^ IGM) as the most valid approximation to neutral 
fraction. At WHIM temperatures, this closely follows the 
CIE neutral fraction, but diverges quickly at T < 10^ K. 
The model ionization parameter U = n^/nn = 0.01 is 
typical of the \ow-z IGM, but large changes in the model 
ionizing field will produce only small changes in fm{T) 
at WHIM temperatures. Figure 1 shows model curves 
for ,fm{T, U) based on photo-thermal CLOUDY models 
with log [/ = -2 and log f/ = -1: thev differ by - 0.3 
dex at T 10^ K, but by only ~ 0.1 dex at T > lO^-^ K. 

Path length AX{Nui) for each sight line is calculated 
as described in Section 3.1. Of the 63 Probable and 
Possible BLA candidates, most are strong enough that 
AX{Nm) ~ AXjnax- However the survey of the weaker 
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Fig. 8. — Distribution of input and corrected 6- values. Input sam- 
ple (dotted) is corrected based on randomly-selected bcOG/''LW 
ratios from DS08. Solid line and shaded boxes show the median 
and ±lcr distribution of lO** Monte-Carlo simulations per absorber. 
Any BLA candidate corrected to 6 < 40 km is dropped from 
the calculation of SIbla performed later in this paper. 



BLAs is only ^ 80% complete. We correct for complete- 
ness by dividing column density iVni by the correspond- 
ing completeness in their respective data (all correction 
factors were between 0.8 and 1.0). The BLA mass frac- 
tion is calculated by modifying equation 6 to 



f^BLA = 



CPcrit 



• (8) 



Given the uncertain relationship between 6obs and 6t 
for any given absorber, we apply a statistical procedure 
to better ascertain thermal line widths. First, we assume 
that 6cOG = '^T where available. For the remainder, we 
perform a Monte-Carlo simulation, correcting each ob- 
served Probable and Possible BLA linewidth as follows. 
Each candidate is "corrected" using a randomly selected 
bcoc/biM ratio from the 138 absorbers in DS08 with 
well-determined 5coG (Fig. 2). Each BLA candidate is 
simulated 10** times and the resulting distribution of hcoi 
is shown in shown in Fig. 8. The survival fraction of 
an individual BLA candidate (i.e., 6cor > 40 km s~*) 
is (62 ±5)% although this includes seven &coG > 40 
km s^^ absorbers which are not corrected. Roughly 50% 
of the non-COG BLA candidates survive the correction 
process. 

The corrected linewidths are used to calculate Orla for 
each Monte-Carlo simulation. The gas temperature is as- 
sumed to be T = 60 fe^oj. K {b in km s^*). Neutral fraction 
as a function of temperature fm{T) is determined from 
a CLOUDY simulation as discussed above, and the total 
hydrogen column density is calculated Nu = Nm/fm for 
each absorber. Summing over all BLAs, we derive flnLA 
according to equation (8). The full Obla distribution for 
10** simulations (Fig. 9) shows a median and ±la value 
of of Qbla = 6.3t X 10-3 or a baryon fraction of 
^BLA/^b = Ull%. Varying some of these assumptions 
changes JIbla as discussed below. 

4.1.1. Systematic Uncertainties and Trends in Qbla 



15 



0^/0,^%) 




30 



logU = -2 
logU = -l 
RFB06 param. 

R04 

R06: secure BLA 
R06: total sample 
L07: CIE 

L07: RFB06 param. 



0.0025 0.0035 0.0050 



0.0070 0.0100 



0.0150 



0.0250 



Fig. 9. — Distribution of Obla values from 10^ Monte-Ccirlo 
simulations (histogram) with median value (filled star) and itlcr 
range (horizontal line) fJeLA = S.Sj^g g x 10~^ ^ro- Various other 
assumptions are also shown for comparison. Using neutral frac- 
tions based on a CLOUDY model with log U = —1 (rather than 
log U = —2) shifts the simulated Hbla distribution to larger values 
f^BLA = (9.6lj'J) X 10^''' (open star). Using the Richter, Fang, & 
Bryan (2006) /hi paramctrization yields SIbla = {^^■S't\'^)x 10~^ 
(open square). Assuming 6l\v = f'T gives a solution ~ 2<t higher 
than the simulations (open diamond). The Richter, Fang, & Bryan 
(2006) assumption that 6t = 0.9 X 6lw (open diamond) gives 
nearly the same value we get from our statistical treatment. Lit- 
erature Hbla values and upper limits from Richter et al. (2004), 
Richter et al. (2006), and L07 are shown for comparison. The 
equivalent baryon fraction, flBi^A/^by is plotted on the top axis 
for comparison. 



There are a number of unccTtaintic^s and assumptions 
present in our JIbla determination whieh ean affect 
f^BLA in several ways. Since we are inferring a total 
hydrogen column density from an observed neutral trace 
component, both 6- values and A/hi have a very large lever 
arm with which to act on the total baryon count. 

First, we assume that the COG-determined (or sta- 
tistically corrected) line width is entirely due to ther- 
mal broadening. As discussed above, even a noise-free 
COG can overestimate a true line width by 20 — 50% 
in the case of undetected, blended components. If we 
assume all &cog values overestimate 6t by 20% (as 
discussed in Section 2.1), the post-correction survival 
fraction of BLAs becomes much lower (34 ± 4%) and 
f^BLA = (2.4+Q 4) X 10~^, more than a factor of two lower 
than our assumed 6coG = &t value. 

Second, since IGM absorbers are thought to be quite 
large in extent and generally unvirialized, there is po- 
tentially a Hubble expansion between one side of the 
cloud and another. Using a typical absorber scale of 
350 kpc (Danforth et al. 2006; Danforth & Shull 2008), 
we might expect a differential Av ~ 24/i7o km s^^ be- 
tween the; two sidcis. This would add in quadrature with 
the other non-thermal line-broadening effects. Running 
the simulation with this additional non-thermal correc- 
tion in place yields a lower mass fraction from BLAs: 
^^BLA = 3.7^1;° X 10-3. 

Third, it is unknown how many of our perceived broad 



lines have a narrow component to them, and conversely, 
how many blended Lya forest lines contain a weak, broad 
component dominated by a narrow absorber. The BLA 
column density is overestimated for the former lines and 
uncounted for the latter. Future observations with high 
S/N (> 30) from the Cosmic Origins Spectrograph may 
disentangle some of these lines, but many blended, mul- 
tiphase systems will likely remain permanently entangled 
(see Fig. 3), leaving the final BLA baryon census uncer- 
tain. 

We assume that the detectability of a line of a partic- 
ular width/depth is a function purely of the local signal- 
to-noisc ratio of the data. Based on this assumption, 
the data are not less than 80% complete for the range 
of BLA candidates. If this assumption is optimistic and 
the data is actually incomplete to ~ 50% in some cases, 
the numerator in Eq. 8 may rise by as much as ^ 60% 
for some absorbers in the sample. The contribution to 
^^BLA from these terms would increase accordingly, but 
we expect the correction to the total Obla sum to be 
minor. 

In this work, we have considered only absorbers with 
^LW ^ 40 km s^^. However, the distribution of 
bhw/bcoG values (Figure 2) shows that 12% (20/164 
COG solutions) have 6lw/^cog < 1- Correcting narrow 
lines by ratios less than unity will result in a statistical 
line broadening and could, in principle, create additional 
BLA candidates. A closer look at the data suggests that 
"scatter-up" is a small effect. The smallest 6-ratio in 
the DS08 distribution is 6lw/&cog ~ 0.7, so only ab- 
sorbers with 28 < 6lw < 40 km would be broadened 



to 6cor > 40 km 



There are 82 Lyct absorbers in 



DS08 in the seven sight lines examined here. On aver- 
age, a maximum of 82 x 0.12 « 10 absorbers would be 
scattered up. However, due to the opposing trends in 
6lw and &lw/^cog distributions, the likely number is 
smaller (4 — 5). Additionally, we have not scrutinized 
the full sample of 30 < 6lw < 40 km s^^ absorbers for 
obvious multiple component structure, which can only 
reduce the likely number of additional BLAs. Thus, we 
expect these "scattered-up" narrow lines to be a small 
correction to Orla- 

4.2. Comparison to Previous Work 

Previous BLA studies have approached the issues of 
line identification, line width, and /hi in different ways 
and most have used some or all of the same datasets 
studied here. It is instructive to look at the assumptions 
in and results from these studies to see how different 
systematics affect dAf /dz and Obla- All of these values 
are shown in Fig. 9 in comparison to our results. 

Richter et al. (2004) study - 9 BLAs in the 
PC 1259-^593 sight hue. They find (cW /dz)BhA « 23 
and calculate /hi in CIE from inferred temperature. 
They note that 6lw overestimates 6t and thus quote only 
an upper limit Obla < 3.3 x 10~^h^Q, roughly half of 
our result. However, this is a single sight line and the 
difference may be explainable by cosmic variance. 

Richter et al. (2006) analyze H 1821+643 and 
PC 0953-1-415 and bring in the results from PG 1259-^593 
(Richter et al. 2004) and PG 1116-^215 (Sembach et al. 
2004). They find 20 and 49 BLAs in their "se- 
cure" and total samples, with {dJ\f/dz)BLA = 22 and 
{dM /dz)BLA = 53, respectively. Assuming CIE, they 
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Fig. 10. — Examples of four Possible BLAs toward PKS 0405— 123 measured in STIS/E140M data which can either be confirmed or 
refuted based on S/N 50 COS/G130M Early Release Observations. The bottom spectra in each panel show the STIS/E140M data used 
in our analysis to derive consensus and b values (listed). The upper spectra are COS data with prominent lines labeled and offset upward 
by 1.5. Panel (a) shows a BLA candidate which has essentially the same profile in both datasets. Panels (b) and (c) show absorbers which 
show a significantly different structure in the high-S/N COS data. Panel (d) shows a reported Possible BLA which does not appear in 
the COS spectrum. All data are smoothed to their respective resolution elements. The consensus line fit of the STIS data is marked by a 
dashed line in both profiles. 



find Obla > 2.7 x lO'^ (3.8 x IQ-^) hj^, but the value 
becomes unphysically large (> fib) if CIE + photoion- 
ization is assumed (dashed curve in Figure 1). 

Lchner et al. (2007) compile data from seven sight 
lines including four from previous work (PG 1259+593 
(Richter et al. 2004), HE 0226-4110 (Lehner et al. 2006), 
HS 0624+6907 (Aracil et al. 2006b) and PG 1116+215 
(Sembach et al. 2004)), two from as- yet-unpublished pri- 
vate communication (PG 0953+415, Tripp et al. in prep, 
and H 1821+643, Sembach et al. in prep), a detailed re- 
analysis of PKS 0405— 123. They perform quality and 
column-density cuts on 341 Lya lines and end up with 
- 60 BLAs for a density (dTV /dz)BLA = 30 + 4. To 
deal with blended components and non-thermal broad- 
ening, L07 randomly eliminate 1/3 of their BLAs with 



40 < 6 < 65 km (all absorbers with 6 > 65 km s~^ 
are kept). Nonthermal broadening is assumed to be 
^ 10% of the total line width (Richter, Fang, & Bryan 
2006), so temperatures are calculated assuming fex ~ 
0.9feobs- They calculate JIbla assuming both CIE and 
Richter 's CIE + photoionization parametrization and 
find SIbla = 3.6 x 10~^ /ifg^ and 9.1 x 10^^ h^^, respec- 
tively or ^BLA/^b = 8% and 20%, for each case. 

At higher redshift, Prause et al. (2007) observed BLAs 
in the spectra of five 1.34 < z < 1.94 AGN using 
both the near-UV STIS/E230M grating on HST and 
the ground-based UV Echelle Spectrograph at the ESQ 
Very Large Telescope. In total, they found 9 good BLA 
candidates and an additional 29 tentative cases in the 
redshift range 0.9 < z < 1.9. They derive a value of 
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^^BLA = (2.2 ± 0.1) X 10^3 for the 9 good BLA can- 
didates and Obla = (14 ± 2) x 10~^h^Q for the entire 
sample. 

4.2.1. Broad Absorbers in Cosmic Origins 
Spectrograph Observations 

Late in the analysis process, we obtained observations 
of PKS 0405-123 by the HST/Cosmic Origins Spectro- 
graph (Green et al. 2010; Osterman et al. 2010) as part 
of its public Early Release Observations (ERO) pro- 
gram. These data were obtained in the G130M grat- 
ing (1132 A < A < 1468 A) with a nominal resolution 
R ~ 20,000 (At; = 15 km s'^). Seven orbits in to- 
tal (~ 17 ksec) were devoted to ERO observations, with 
three taken prior to when an accurate focal alignment 
was achieved and four afterward. A close examination 
of the data shows no difference in line profiles for any of 
the ISM or IGM lines of interest, so observations from 
all seven orbits were aligned and coadded. 

The resulting spectrum is of exquisite quality with 
S/N > 40 per nominal seven-pixel resolution element 
at most locations. Owing to the differences in detec- 
tor technology and the different grating positions used 
in the observations, the COS data are free of much of 
the fixed-pattern noise that plagues the corresponding 
STIS/E140M observations. This, coupled with the very 
high S/N. makes COS ideal for verifying the BLA can- 
didates discussed in this paper. 

Eighteen of the BLA candidates toward PKS 0405-123 
measured in STIS/E140M data are also covered in the 
COS observations. Five of these cases are blended com- 
ponents of strong absorption systems, which are difficult 
to confirm or refute, but the COS data is consistent with 
the STIS observations. Of the 13 weaker absorption lines, 
seven show a good match between STIS and COS data. 
However the STIS-measured line profile is substantially 
different in four cases and missing altogether in another 
three. Figure 10 shows several examples where BLA can- 
didates can be either confirmed or refuted based on COS 
observations. 

This is a good demonstration of the importance of 
fully understanding the instrumental effects of a par- 
ticular spectrograph. While the COS data are not free 
of fixed pattern noise, it is independent of that in the 
STIS observations. Furthermore, it is evidence that the 
~ lOx sensitivity increase of COS over STIS will revo- 
lutionize the study broad Lya absorbers, both through 
higher S/N and sensitivity to numerous, fainter targets. 
Since this work relies on a uniform analysis of consistent 
datasets by independent groups, we do not change any 
of our BLA designations in light of new, higher-quality 
data from COS. However, we note the results of our COS 
absorber verification where possible in the individual ab- 
sorber comments in the Appendix. 

5. CONCLUSIONS AND SUMMARY 

Broad Lya absorbers are a potentially powerful 
method of measiiring the extent and distribution of gas 
at T > 10'^ K in the intergalactic medium without re- 
lying on metal enrichment. The small hydrogen neutral 
fraction even at WHIM temperatures will result in broad, 
shallow H I profiles that can be translated into temper- 
atures and total hydrogen column densities. Unfortu- 
nately, BLAs present some observational challenges and 



ambiguities. Broad, shallow absorbers are difRcult to de- 
tect in data of only moderate S /N and in multi-phase sys- 
tems. Detected BLAs are strongly biased toward cooler 
temperatures where lines are relatively narrower and the 
neutral fraction is higher. Furthermore, identification of 
bona fide BLAs relies crucially on the line component 
structure, knowledge of instrumental features and an ac- 
curate continuum definition. 

We attempt to work around these problems by inde- 
pendent reduction and analysis of seven AGN sight lines 
containing 119 purported broad Lya lines reported in the 
literature (mainly DS08 and L07 and sources therein). 
We assign consensus values for column density and line 
width based on two independent analyses by two dif- 
ferent research groups. The purported BLAs are split 
into three qualitative categories based on the consen- 
sus lincwidth, detailed analysis of the absorption pro- 
file, and other factors. Probable BLAs (15 systems) are 
those showing 6lw > 60 km s~^ and no obvious asym- 
metry or component structure. Systems with curve-of- 
growth determined linewidths 6coG > 40 km s~^ are 
also deemed Probable BLAs. Possible BLAs (48 sys- 
tems) are those with 40 < 6lw < 60 km s^^, those ab- 
sorbers with potential component structure or asymme- 
tries, or some plausible reason to doubt their identity as 
BLAs. The remaining systems (56) are ruled out as being 
BLAs for a number of reasons: alternate line identifica- 
tion, 6lw < 40 km s^^, probable component structure, 
or simply failing to appear in our reduction of the data. 

Taking all of the probable BLAs and 50 ± 50% of the 
possible category, we see 39 ± 24 BLAs along a total red- 
shift pathlength Az = 2.193 (AXtot = 1-773) in the 
seven AGN sight lines surveyed. This gives a detec- 
tion frequency of dJ\f /dzsLA = 18 ± 11 {dM /c?Xbla = 
22 ± 14). This frequency is similar to that of O VI, an- 
other potential WHIM tracer with {dAf /dz)ovi = 15 ± 3 
(Danforth & ShuU 2008; Tripp et al. 2008), though the 
BLA frequency has considerably greater uncertainty. In- 
deed, while the detection or non-detection of highly ion- 
ized ions (O VI, N V, Ne VIII, etc) was not taken into 
account in our BLA categorization, 40% of the probable 
BLAs and 20% of the combined probable-plus-possible 
samples show reasonable O VI detections. The incidence 
of O VI detections in narrow Lya lines is ~ 15%. 

The relationship of WHIM to galaxies is another key 
area of interest. It is likely that small galaxies with weak 
gravitational fields are important for IGM heating and 
enrichment (e.g., Stocke et al. 2004). Unfortunately, sur- 
veys for low-luminosity galaxies tend to be unreliable at 
redshifts greater than a few hundredths. However, we 
found eight BLAs in regions of fairly complete galaxy 
surveys (L > 0.2 L*) that provide nearest-galaxy dis- 
tances of 60 kpc out to nearly 3 Mpc. The detection 
of O VI in conjunction with broad Lya was highly corre- 
lated with galaxy distance as the four BLAs with d < 600 
kpc showed O VI detections at some level, while the four 
BLAs at > 600 kpc appear to be free of O VI ab- 
sorption. This result gives significant support for BLAs 
probing gas that O VI surveys do not detect. Using the 
above O VI detection statistics, our BLAs suggest that 
~ 80% of the baryons are not already accounted for in 
O VI surveys. 

A main and crucial uncertainty in BLA surveys is dis- 
entangling thermal and non-thermal line broadening, as 



18 



measured line width overpredicts the thermal linewidth, 
often by an unknown amount. Previous studies have rec- 
ognized this phenomenon and dealt with it in a variety of 
ways, including arbitrarily throwing out some fraction of 
measured BLAs from a sample or scaling measured line 
widths by some uniform factor. We approach the prob- 
lem by looking at curvc-of-growth 6- values (ftcoc): the 
doppler b-parameter measured from a single line (typ- 
ically Lya) overpredicts that from a full COG analysis 
(Danforth et al. 2006) by a factor of ~ 1.5 ± 0.9 (though 
we also show how even a COG can overpredict the true 
linewidth). Since most broad Lya systems arc fairly 
weak, confirmation in higher-order Lyman lines is usu- 
ally not possible. Instead, we statistically correct for the 
single-line 6-value overprediction based on the observed 
&Lw/fecoG distribution in DS08. We find in our simu- 
lations that 55 ± 5% of the reported BLAs survive the 
correction process with fecor > 40 km s^^, corresponding 
to T > 10-^ K. 

From temperatures derived from line widths, we can 
estimate the neutral fraction of a particular absorber 
and hence the total hydrogen column associated with 
a particular H I detection. The total amount of gas 
at T > 10^ K traced by these BLAs can then be esti- 
mated as a fraction of the closure density of the Uni- 
verse. We use a Monte-Carlo simulation to statistically 
correct the observed BLA linewidths. Our median value 

of riBLA = 6.3tJ-i X 10-3 /l-l (f^BLA/^^fa = 14^^%) 

based on 10^ Monte-Carlo simulations of each BLA. Since 
~ 80% of these BLA baryons are not in O VI systems, 
the combination of O VI -I- BLA WHIM searches can ac- 
count for r2vvHiM/f^5 ^ 20%. This includes ^ 8% from 
O VI (DS08) and an additional ~ 12% in metal-poor 
BLAs. It is clear that systematic uncertainties involved 
in BLA surveys are comparable to or larger than the 
statistical fluctuations from cosmic variance among the 
sight lines, and more work must be done to understand 
the individual systems. 

To illustrate the importance of methodology and in- 
dividual systems to the baryon census, we compare our 
result to that of L07. L07 used the same set of AGN 
sight lines used here and a more inclusive set of BLA 
candidates. They assume 6t = 0.9&lw as in Richtcr, 
Fang, & Bryan (2006), calculated hydrogen neutral frac- 
tion via a CIE assumption (very similar to the CLOUDY 
log N = —2 model used here), and randomly eliminated 
1/3 of the 40 < 6 < 65 km s"'^ absorbers which gave 
a result of Obla = 3.6 x 10"^ h^Q or ^BhA/^b = 8% 
or 60% of our value. If, instead, /hi is based on the 
CIE -|- photoionization model of Richter, Fang, & Bryan 



(2006), the L07 value rises to f^BLA = 9.1 x lO"'"' hj^^, or 
~ 40% larger than our result (see Fig 9). 

Future observations of low-redshift BLA systems with 
the Hubble Space Telescope/ Cosmic Origins Spectro- 
graph (COS) will improve the BLA census in several 
important areas. First, five of the seven sight lines stud- 
ied here (HE 0226-4110, PG 1116-1-215, PKS 0405-123, 
PG 0953+415, and PG 1259+593), as well as over a dozen 
other AGN sight hncs observed with STIS/E140M will 
be observed by COS as part of the Guaranteed Time 
Observations. Some of the scheduled GTO observations 
are high-S/N spectra of BLLac objects chosen specifi- 
cally to search for BLAs against non-thermal power-law 
continua. More AGN observations are scheduled in sev- 
eral large HST Cycle 17 Guest Investigator programs 
(Pis: Tripp, Tumlinson). Observations of the same tar- 
gets by different instruments will help sort out real BLAs 
from instrumental features. Furthermore, the exquisite 
S/N expected from COS data, both for previously ob- 
served and new sight lines, will be crucial in determin- 
ing line profiles, identifying blended systems, and finding 
the expected population of weak, broad systems. Late in 
our analysis process, we obtained high-S/N COS/G130M 
observations of PKS 0405— 123 which shed considerable 
light on individual BLA candidates. We discuss this fur- 
ther and show several examples of possible BLAs ob- 
served with much higher S/N in the Appendix. While 
some BLA Candidates are confirmed by the COS data, 
there are some significant differences which suggest that 
the actual number of BLAss less than catalogued herein. 
Finally, the high sensitivity of COS compared with STIS 
will allow a much larger pathlength of the low-z IGM 
to be surveyed, increasing our statistics on intergalac- 
tic absorbers ranging from metal-line systems to BLAs. 
Increases in both the O VI and BLA catalogs will un- 
doubtedly occur. It will be very interesting to see where 
the WHIM baryon census stands ten years hence. 
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APPENDIX 

NOTES ON INDIVIDUAL ABSORBERS 

Wc present here more detailed comments on individual BLA candidates than is present in Table 2. Each system is 
identified by sight line, redshift, and our BLA classification (A=Probable BLA, B=Possible BLA, C=not a BLA). We 
add a fourth category here (X) for narrow Lya systems which are adjacent to or blended with BLA candidates. See 
Table 2 

HE 0223-4110 

z = 0.06083 (B):— Strong line detected in Lya,/3 (M^Lya = 557 ± 20 mA, Wby/j = 313 ± 113 mA) with suggestion 
of broad wing on blue side and possible component structure. No metal lines to clarify component structure. COG 
for full system gives 6coG ~ 36 km s~^ with large uncertainties. We list this as a possible BLA since there the broad 
wing on the blue side suggests a broad component. 

z = 0.16339 (B): — Very strong Lya system. Ly/3 is blended with Galactic Si II absorption, but appears narrow. 
Ly7 absorption is present but weak. COG solution doesn't converge and there are signs of multiple components. No 
metal lines to illustrate component structure. 

z = 0.20700 (B): — Saturated, broad Lya profile with possible weak, broad component seen in the wings. Ly/3,7 
show asymetric profiles. Strong O VI, C III, Si III (DS08) and Ne VIII (Savage et al. 2005) detections as well as possible 
S IV. We adopt the Savage et al. (2005) measurements for the BLA component, but list the system only as a possible 
BLA since the component structure is ambiguous. 

z = 0.30930 (B): Strong, broad system with probable component structure in Lya,P,j: WLya = 388 ± 7 mA, 
WLy/3 = 130 ± 3 mA, WLy-y = 33 ± 3 mA. Multi-valued COG solution: bcoG = 34 ± 2 km s'^ {a, (3) or 6coG = 
44 ± 6 km s~^ (a, 7). Due to this ambiguity and the likely component structure, we list this as a possible BLA. No 
metal lines. 

z = 0.38420 (C): Strong system with ambiguous component structure in Lya. However, Ly/3 profile shows two 
clear components. Two or three components seen in Ly/3 profile despite poor S/N. COG on stronger, red, component 
gives bcoG ^ 31 km s~^. 

z = 0.39641 (B). z = 0.39890 (B), z = 0.40034 (B): — Possible continuum ripples shown in Figure 5 and discussed in 
the text. In particular, a very similar feature appears in the data at 1702. 5A in the PG 0953+414 data at Zabs > -^agn 
which looks similar to the z = 0.40034 feature toward HE 0223—4110, however other spectra at higher S/N do not 
show this feature. 

PKS 0405- 123 

z = 0.03196 (B):— Broad Lya profile (W^Lyo = 108 ± 22 mA, 6lw = 57 ± 8 km s^^), but weak, narrow Ly/3 
(M^Ly/3 = 50 ± 6 mA, 6lw = 35 ± 4 km s~^). COG is poorly constrained: 6cOG > 10 km s~^. Possible BLA based on 
Lya alone. COS observations are consistent with the STIS/E140M data. 
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z = 0.05896 (C): — Noisy data suggests two weak, narrow components each with 6 ~ 20 km s ^. COS spectrum 
confirms the two components with 6lw = 23 and 27 km are not BLAs. 

z = 0.07218 (C):— Very marginal feature which is not clearly Lya. Identified as C III 2; = 0.3340 by DS08. 
Absorption feature confirmed in COS data as real with &lw = 30 ± 3 km s^-'^. 

z = 0.07523 (B):— Very weak detection measured by both W06 (6lw = 56 km s'^) and L07 (6lw = 48±20 km s'^). 
High-S/N COS observations show no absorption at this position iyV < 11 mA). 

z = 0.08139 (A ):— Moderately strong system with 6lw ~ 53 km s'S WLya = 261±16 mA. Ly/3 shows W = 90±8 mA 
but is blended with Lye z = 0.18269 (WLya ~ 670 mA). If this entire absorption feature is taken as Ly/3 z = 0.08139, 
then 6coG > 22 km s~^. However, if only part of absorption is taken as Ly/3 {Whyfs ~ 45± 10 mA), the COG linewidth 
is large but uncertain: 6coG = 74 ± 30 km s~^. Given this ambiguity and the broad, clean Lya profile, we list the 
absorber as a probable BLA. 

z = 0.09659 (B): — Strong Lya system (Whya ~ 500 mA) with excess in blue wing indicating possible broad 
component. DS08 find 6coG = 36±4 km and W06 find ^lv^^ = 37± 1 km for the entire system. However, L07 
fit a broad component to the system and we take this as a possible BLA due to the ambiguous nature of the fit. DS08 

report O VI detected in both lines of the doublet associated with the system, though it is unclear to which component 
it may correspond. COS observations show consistent line profile, but BLA component cannot be verified. 

z = 0.10298 (C): — Very marginal feature which is not obviously a single absorber. W06 measure 6lw = 60 ± 
14 km s~^. COS data shows a clear blend of two absorbers 6lw = 27 ± 5 and 42 ± 3 km s~^, thus no BLA is confirmed 
in this system.. 

z = 0.10419 (C): — Very marginal detection not confirmed by DS08. COS observations show no absorption at this 
position {W < 10 mA). 

z = 0.13102 (B): — Moderate absorption feature with broad wing on red edge suggesting possible BLA component. 
COS observations show absorption at this wavelength (6lw = 35 ± 1 km s"""^, log N = 13.18 ± 0.01), significantly 

narrower than that observed by STIS. 

z = 0.13377 (B): — Noisy data with moderately broad absorption. COS observations show a clear absorber at this 
location (6lw = 53 ± 2 km s~^, log N = 13.10 ± 0.02) but significantly shallower and broader than the STIS profile. 

z = 0.13646 (C): — Marginal line with likely components. Single-component fit from W06 gives 6lw = 49±8 km s~^. 
COS confirms double component structure, 6lw = 23 ± 4, 26 ± 2 km s~^. 

z = 0.13924 (C): — No detection in L07. Marginal feature in DS08 and W06. Not clearly a single component even 
if real. COS observations show two absorbers offset in velocity, but nothing corresponding to the reported line profile. 

z = 0.15304 (C): — Strong, broad line (WLya = 263 ± 10 mA) identified as Lya by L07, but corresponding Ly/3 
nondetection (WLy/3 < 24 mA) is inconsistent with a single Lya absorber. Unknown line identification. COS confirms 
line profile. 

z = 0.16121 (C): — Strong absorber with several clear components. Strong metal lines likely make up the fianking 
components leaving little room for a broad Lya core. COS confirms general profile. Thus COS does not confirm this 
BLA. 

z = 0.16678 (B) and z ~ 0.16714 ( C): - Strong Lya complex (T4^Lya ^ 650 mA) with an asymmetry on the blue wing. 
Line can be fit with broad a broad and narrow component, but decomposition is suspect. Ly/3 {Wi^yp ^ 450 mA) shows 
no sign of the broad, blue component. Strong system (z = 0.16714) shows obvious component structure in various 
metal lines. The system shows strong O VI and N V absorption (DS08), but it is unclear with which component it 
is associated. COS observations show a consistent line profile, but cannot confirm or deny the presence of a BLA 
component. 

z = 0.17876 (B): Asymmetric Lya line with very marginal Ly/3 detection. COG analysis by DS08 gives &coG = 
18^7^ km s~^. Apparent Optical Depth (AOD) line fits to Lya line gives much broader profile (6lw ~ 60 km s~^; W06 
measures 6lw = 58 ± 6 km s~^). Well fit by an offset broad+narrow pair (6lw = 57, 14 km s~^), but inconclusive. 
We take the AOD measurements as concensus in this case. COS observations confirm the general profile shape of this 

absorber and consistent line measurements. 

z = 0.18269 (B): — Strong H I system (WLya ~ 690 mA) with at least two components, possibly more. Ly/3 is 
blended with Galactic Lya. Ly7 shows strong, moderately asymmetric profile iW-Ly-y = 139 ± 11 mA). COG gives 
''COG = 49^7^^ km s^^, but probable component structure relegates this to a possible BLA. System shows several 
O VI components. COS observations confirm the line profile of this strong absorber, but cannot confirm or deny the 
presence of a BLA. 
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z = 0.19086 (B): — Weak, noisy absorber with possible component structure. COS observations do not show an 
absorption line at this position {Wr < 22 mA). 
The COS ERO spectral coverage (1132 < A < 1468) ends at this maximum redshift. 

z = 0.24513 (B) and z = 0.24553 (X): — Pair of well-separated Lya absorbers. System at z = 0.24513 appears 
narrow (6lw = 30 ± 5 km s~^) but both L07 and W06 measure 6 « 55 km s~^. We adopt a consensus b and N 
intermediate between the two measurements to account for possible continuum and reduction differences. 

z = 0.32500 ( C): Marginal detection, however there is a correlation between features in Lya,/?. Probably multiple 
components. DS08 report 6coG = 21 km s^^, however we find this to be more reasonable as a lower limit and the 
COG solution doesn't fit the absorption profile very well. W06 report 6lw = 81 ± 11 km s~-^. Higher S/N data on 
this absorber should prove interesting. Given the marginal nature of this feature and the likely component structure, 
we list it as a non-BLA. 

z = 0.35092 (B), z = 0.35149 (X):— Strong Lya system flanked by Si III z = 0.36079 and Lya z = 0.35149. Lya,/3 
cinvc of growth gives 6coG = 56 ± 15 km s^^ (DS08), but Lya, 7 curve of growth gives 6coG = 27 ± 5 and seems 
to better match the profile. W06 and L07 measure bLya = 38 ± 2 km s~^ and ^Lya = 40 ± 8 km s~^, respectively. 
Both Lya and Ly/3 show asymmetric profiles suggesting multiple components. We accept the Lya measurements as 
consensus in this case and list the system as a possible BLA. 

z = 0.36150 (C) and z = 0.36079 (X): Strong, blended Lya systems with an ambiguous component structure 
and noisy data. Two-component fit gives &lw = 75 ± 10, 39 ± 5 km s~^ for the z = 0.36150 and z = 0.36079 
components, respectively. The z = 0.36150 component, however, doesn't appear in Ly/3 (VFLy/3 < 9 mA); based 
on the Lya component, the expected Ly/3 line should be WLy/3 ~ 50 mA. Meanwhile, a Ly/3,7 curve of growth fit 
to the the z = 0.36079 component (ignoring the noisy, possibly blended Lya profile) gives 6cOG = 25 ± 3 km s~^, 
log N = 15.18 ± 0.13 cm~^. We conclude that neither system is a BLA. 

z = 0.40886 (B):— Strong Lya system (V^Lya = 430± 10 mA) at the long-wavelength, noisy end of the STIS/E140M 
range. Ly/3,7 also present (W^Ly/3 = 126 ± 15 mA, W^yj = 68 ± 10 mA), both with asymmetric profiles. Lya,^ and 
Lya, 7 COG solutions inconsistent but generally 6coG ~ 40 km s~^. Possible BLA, but probable component structure. 

HS 0624+6907 

z = 0.04116 (C): — Weak feature consistent with Si III z = 0.06346, a system in which many other ionic lines are 
also detected. 

z = 0.05437 (B), z = 0.05515 (B), and z = 0.05484 (C):— Strong Lya system (WLya ^ 450 mA) with several 
plausible BLA subcomponents. The z = 0.05437 line appears as a very marginal, noisy absorption on the blue wing 
while the z ~ 0.05515 component is an excess on the red side of the main absorption line. The strong central component 
at z = 0.05484 is fit well by a 6lw ^ 35 km s^^ component despite the &cOG = 45lg^ km s^^ value reported in DS08. 
The narrow feature at 1283.15 A can be unambiguously identified as Si III z = 0.0635. The BLA nature of these lines 
relies on the details of the component structure which is quite ambiguous and difficult to disentangle. 

z — 0.06346 (C): — Strong Lya system with measured &lw > 40 km s~^ but 6cOG < 40 km s~-^ with O VI detected 
in both lines of the doublet (DS08). H I profile looks symetric, but there are two clear components in Si IV, C IV, 
Si III, and S III. 

z = 0.13597 (B):— Weak double profile. Feature at 1381.4A is Ly/3 z = 0.3468 but broader feature at 1381.0A is a 

plausible BLA with low-significance O VI detections at both 1032 and 1038A. Possible BLA based on the measured 
line though the large error bars and O VI detection coiild argue for this being a Probable BLA. 

z = 0.30994 (B): — Broad, weak absorption feature identified as Lya by both L07 and DS08. However considerable 
continuum uncertainty yields highly incompatible b, N solutions. Possible, low-significance O VI detected as well. In 
any case, profile is asymmetrical and not obviously a single component. 

z = 0.31045 (C), z = 0.31088 (C): — Broad, asymmetric features which are not obviously real. Possible continuum 

differences. 

z = 0.31790 (B): — Moderate, narrow Lya line with asymmetry on the red wing and good O VI detections in both 
lines of the doublet. Possible multiphase system and/or continuum fit uncertainties. Possible BLA. 

z = 0.32089 (C): — Strong, moderately-narrow line is a poor fit to DS08 COG solution (6coG = ^4+11 km s'^). 
Lya absorption is anomalosly large given observed Ly/3,7, however there are no other obvious contributors to the Lya 
line strength. We confirm the &lw = 31 ± 1 km s~^ solution of L07 for this absorber. 

z = 0.33976 (A): — Strong, broad system with no obvious components and detections in Lya,/3,7. Good COG 
solution gives 6coG = 41 ± 6 km s~^. O VI detected in both lines of the doublet (DS08). 
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PG 0953+415 

z = 0.04382 (C):— Reported as O VI A1032 at z = 0.22974 by Tripp et al. (2008). Interestingly, this is a "blind" 
O VI system with little or no accompanying Lya absorption (WLya < 7 mA). 

z = 0.05879 (B): — Narrow line with a possible excess in line wings suggesting a multiphase system. Very ambiguous 
decomposition. Wc take the L07 values as the consensus. 

z = 0.12784 (C): — Claimed broad blue wing on a narrower Lya system. We see no good evidence for this system 
in our reduction of the data. 

z — 0.17985 (B): — Weak single or double absorption. Potentially fit as single BLA or two narrow components 
(6lw ~ 20). 

z = 0.19126 (B): — Asymmetric Lya profile with uncertain decomposition. Formal AOD linewidth is 6lw ~ 
66 km s^^ however, two components (6lw = 40, 32 km s~^) provides a better fit. We adopt 6lw = 40 : km s~^ and 
log N = 13.3 cm^^ as consensus values. 

z = 0.19361 (C): — Strong narrow Lya system with very weak Ly/3,7 counterparts. DS08 report a Lya,/3 solution 
(^COG = 24^3 km s~^, log N = 14.15 ± 0.04 cm~^) which is stronger than the data can easily support. We measure 
6lw = 37± 2 km s~^ instead and adopt that value here. 

PG 1116+215 

z = 0.08587 (B), z = 0.08632 — Pair of very shallow, ripples in the continuum with 6lw ^ 55 and 6lw ^ 34 km s~^, 
respectively, both with considerable uncertainty. Broader feature is possible BLA. 

z = 0.09279 (A): — Very weak, broad system blended with two weak, narrow features (possible Galactic C I). While 
profile looks questionable, it is more obvious when looking at a broader wavelength range. Eliminating the two sharp 
features gives a very plausible BLA. 

z = 0.13370 (A):— BLA with a low-significance O VI A1032 detection (log A^ovi ~ 13.2 ±0.3 cm'^). Potentially fit 
as two components (A ~ 1378.3, 1378.7) but 6lw ~ 50 km s~^ even in this case. 

PG 1259+593 

z = 0.00229 (A): — Lya line appears on the edge of the broad Galactic Lya trough and Ly/J is blended with an O I 
airglow line in FUSE data. Richtcr ct al. (2004) used night-only data to get a good Ly/3 measurement and constrains 
^COG = 44^4 km s^^ which wc adopt here. 

z = 0.04606 (C):— Very strong system with obvious components in O VI, C IV, Si III, Si IV. Richter et al. (2004) 
report two H I components, both with 6coG < 40. 

z = 0.14852 (A): — Strong nearly triangular line profile with detections in Lya,/3. Curve of growth consistent with 
BLA 6coG ~ 57 km s~^ though with large errors. 6Ly/3 ~ 36 km s~^. 

z = 0.19573 (B), z = 0.19620 (X):— Double absorption profile. Shallow, blue feature at 1453.6 A is plausible BLA 
candidate (6lw = 46± 19 km s~^) and was reported as a component to the z = 0.1962 system by Richter et al. (2004). 
Marginal N V detection. 

z = 0.21136 (C) and z = 0.43569— Richter et al. (2004) identify the absorption at 1472.6 A as Ly/3 z = 0.43569 
and high-res FOS spectra confirm Lya line at this redshift. Absorption at 1242.5 A where Ly/3 z = 0.21136 would be 
expected is blended with Lya z = 0.02217 and unrecoverable. Can't confirm. 

H1821+643 

z = 0.11133 (A), z = 0.11167 (X): — Double absorber profile with one broad, one narrow system (6lw = 
52, 31 km s~^). L07 measure this as a single system (6lw = 88 ± 14 km s~^). One BLA likely present. 

z = 0.12147 (A), z = 0.12117 (X): — Strong and highly asymmetric absorber. Red wing Ly/3 is hinted at in the 
data, but at low significance. System shows 6cOG > 40 km s~^, but there are clearly multiple components. Weak 
O VI absorption is reported in the system (Tripp et al. 2001, 2008; Oegerle et al. 2000) aligned with broad Lya wing, 
however, DS08 do not confirm this detection. 

z = 0.14754 (B) and z = 0.14776 (X): — Blended system of two, easily separable components (6lw ~ 42, 20 km s~^). 
Marginal Ly/3 detection for both components but COG solution is very poorly constrained (6coG = 22^55 km s~^) 
for the z = 0.14754 system. 

z = 0.21326 (A):- Strong system with Lya,/3 detections (&cOG = 42^4 km s"""^) though with some hint of component 
structure. Strong O VI detection in both lines of the doublet. The entire absorber may not represent hot gas, but 
there is probably a BLA component. 
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z = 0.22480 (C): — Very strong system (WLya ~ 800 mA) with clear components in higher Lyman hnes. O VI 
detected. Not confirmed by L07, but clearly not a BLA. 

z = 0.22616 (B): — Moderate absorption feature (WLya = 150 ± 8 mA) identified as Lya by L07. However, 
corresponding Ly/3 non-detection shows WLy/3 < 10 mA, inconsistent with a single Lya system. Possible O VI 
absorption (log A^ovi = 13.4 ± 0.1 cm~^) in both lines of the doublet offset ^ +60 km s~^ from the purported H I 
system. Probable multiple components or misidentification. 

z = 0.25814 (B): — Strongly asymmetric profile suggests a multiphase H I system. Profile is reasonably fit by two 
components: 6lw ~ 56, 17 km s~^. Total line shows WLya = 141 ± 3 mA and corresponding weak Ly/3 absorption 
(WLy/3 « 24 mA) gives a curve-of-growth solution 6cOG = 30 ± 7 km s ^ for the combined system. In light of this, the 
component fits seem reasonable and we list this as a possible BLA. 

z = 0.26658 (A): — Strong system with O VI detected in both lines of the doublet. Poorly-constrained COG gives 
^COG = 68 ± 30 km s~^. Individual lines show 6lw = 46 km s~^ and 6lw ~ 40 km s~^. Probable BLA despite 
poorly-constrained COG. 



